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Abstract 

This  study  investigated  methods  of  target  detection  using  Wide-Angle  Synthetic  Aperture 
Radar  (WASAR).  WASAR  uses  multiple  aspect  angle  Synthetic  Aperture  Radar  (SAR)  images  of 
the  same  scene.  The  SAR  images  were  generated  using  a  pre-release  software  package  from  Loral 
Corporation.  The  software  was  able  to  generate  512  by  512  pixel  SAR  images  that  contained  various 
vegetation  returns  which  for  our  purposes  we  classified  as  clutter.  Within  this  clutter,  targets  (M35 
trucks)  could  be  placed  at  random  locations  and  orientations.  The  software  also  had  the  capability 
of  generating  fully-polarimetric  WASAR  images  with  multiple  depression  angles.  This  data  was 
then  processed  and  various  detection  algorithms  tested  to  exploit  the  amount  and  diversity  of 
information  available  from  the  multiple  images.  SAR  images  are  generally  known  to  contain  large 
amounts  of  data  and  WASAR  images  contain  even  more  due  to  the  multiple  images.  Various  pre¬ 
processing  filters  were  analyzed  for  detection  optimization.  These  filters  included:  polarimetric 
averaging,  polarimetric  span,  polarimetric  optimal  weighting,  and  polarimetric  whitening  filter. 
Simple  classical  detection  (thresholding)  algorithms  were  evaluated  using  these  pre-processed  data 
sets.  The  use  of  WASAR  imagery  improved  detection  by  allowing  thresholds  to  be  set  higher  than 
for  simple  SAR  thereby  avoiding  false  alarms  yet  still  allowing  detection  of  the  known  targets. 
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Processing  of  Wide-Angle  Synthetic  Aperture  Radar  Signals  for  Target  Detection 


/.  Introduction 

1.1  Problem  Statement  and  Significance 

This  thesis  investigates  the  detection  of  targets  using  Wide-Angle  Synthetic  Aperture  Radar 
(WASAR).  The  detection  and  identification  of  military  targets  (buildings,  vehicles,  roads,  runways, 
etc.)  is  important  to  the  Air  Force  mission  of  strategic  and  tactical  reconnaissance  and  bombard¬ 
ment.  To  this  end,  recently  available  WASAR  imagery  needs  to  be  exploited  for  detection  purposes. 
Synthetic  Aperture  Radar  (SAR)  produces  high  resolution  imagery  using  a  moving  radar  platform 
to  synthesize  the  resolution  cells.  SAR  has  the  advantage  of  being  relatively  weather  invariant  as 
compared  to  photo  reconnaissance.  As  an  expansion  of  SAR,  WASAR  provides  multiple  viewing 
angles  of  the  same  SAR  image.  Wide-Angle  SAR  improves  on  SAR  by  providing  more  information 
on  a  scene.  By  viewing  a  scene  from  several  different  angles,  it  may  be  possible  to  pick  up  a  “glint” 
at  one  of  those  angles  indicating  the  presence  of  a  target.  Also,  most  modern  SAR  systems  provide 
the  data  in  fully  polarimetric  format  which  will  be  explained  more  fully  later.  Fully  polarimetric 
data  provides  even  more  information  on  a  scene  over  that  provided  by  a  system  capable  of  only 
a  single  polarimetric  return.  The  diverse  information  obtained  from  WASAR  provides  the  target 
detector  with  vast  amounts  of  information  that  should  make  detection  much  easier.  This  thesis 
takes  a  first  look  at  the  use  of  Wide-Angle  Synthetic  Aperture  Radar  imagery  for  target  detection. 

1.2  Background  and  Problem  Areas 

Before  moving  into  the  processing  of  the  WASAR  images,  it  may  be  helpful  to  go  over  some  of 
the  terminology  that  will  be  used  throughout  this  thesis.  Target  detection  versus  target  recognition 
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Figure  1.1  SAR  Image  Creation 

will  be  discussed,  along  with  terms  associated  with  SAR  and  WASAR.  Some  of  the  problems 
involved  in  using  SAR  and  WASAR  will  also  be  discussed. 

The  first  concept  that  needs  to  be  understood  is  Synthetic  Aperture  Radar.  SAR  uses  a 
moving  platform  radar  system  to  synthesize  a  high-resolution  image.  The  resolution  of  the  image 
is  determined  by  the  pulse  width  (range  resolution)  and  the  distance  the  platform  travels  (azimuth 
resolution).  These  two  parameters  form  a  rectangular  resolution  cell  in  space,  as  seen  in  Figure  1.1. 
Through  geometrical  processing,  these  resolution  cells  are  then  combined  to  create  an  image.  Typ¬ 
ical  SAR  imagery  is  shown  in  figure  2.9  on  page  2-15.  This  image  shows  a  forested  region  on  a 
grassy  plain.  Interspersed  within  this  image  are  targets.  The  targets  in  some  instances  look  similar 
to  trees  in  the  image. 

Wide-Angle  Synthetic  Aperture  Radar  (WASAR)  differs  from  normal  SAR  by  generating 
multiple  images  of  the  same  scene.  The  images  are  formed  over  several  target  aspect  angles. 
Usually  these  different  aspect  angles  are  the  result  of  a  very  wide  beam  that  is  divided  into  several 
smaller  beams,  as  shown  in  Figure  1.2.  As  each  of  these  beams  pass  over  the  area  of  interest,  the 
radar  receives  a  first,  a  second,  and  a  third  look  at  the  same  scene  from  different  viewing  angles. 
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The  number  of  looks  depends  on  the  number  ci  narrow  beams  carved  out  of  the  wide  beam.  This 
research  used  seven  aspects.  Currently,  there  are  no  operational  systems  that  use  this  technique; 
however,  the  use  of  WASAR  for  various  reconnaissance  and  surveillance  purposes  is  being  explored 
by  several  agencies  including  the  Air  Force’s  Wright  Laboratory  at  Wright-Patterson  Air  Force 
Base,  Ohio.  One  of  the  problems  with  WASAR  is  the  registering  of  the  multiple  images  with 
respect  to  each  other.  Since  the  images  are  taken  from  different  aspect  angles,  the  coordinate 
system  for  one  image  is  not  the  same  as  that  in  another.  Registering  the  coordinate  systems  is 
important  if  the  images  are  to  be  compared  to  each  other  in  any  way.  Variations  in  the  radar 
platform  can  cause  major  problems  with  this  coordinate  registering.  For  example,  if  the  range  to 
target  distance  varies  from  one  aspect  to  another,  the  ground  plain  resolution  will  not  be  the  same 
from  image  to  image.  In  other  words,  a  pixel  in  one  image  may  represent  2  square  feet;  whereas, 
in  another  it  may  represent  3  square  feet.  This  particular  problem  of  coordinate  registering  can 
pose  quite  a  problem  when  the  radar  system  is  an  airborne  SAR  platform.  In  this  case,  either  the 
aircraft  has  to  fly  a  very  precise  route,  or  the  images  have  to  be  registered  using  some  type  of  signal 
processing  system  that  uses  not  only  the  images  but  flight  data. 
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Additionally,  SAR  and  WASAR  imagery  can  be  generated  using  polarized  data.  Polarization 
is  a  term  used  to  describe  the  electro-magnetic  wave  orientation  of  the  radar  transmit  and  receive 
signals.  Fully  polarimetric  radar  returns  consist  of  three  distinct  modes:  HH,  HV,  and  VV.  The  HH 
term  describes  the  horizontally  polarized  electro-magnetic  (EM)  field  resulting  from  illumination 
of  the  scene  by  a  horizontally  polarized  EM  source.  HV  polarization  describes  the  horizontally 
polarized  EM  field  resulting  from  illumination  by  a  vertically  polarized  EM  source.  In  order  for 
a  radar  data  return  set  to  be  fully  polarimetric,  it  must  contain  all  three  of  these  elements. 
Polarimetric  processing  is  the  combining  of  the  HH,  HV,  and  VV  images  in  some  optimal  fashion 
in  order  to  improve  target  to  clutter  ratios. 

To  this  end,  Holm  [7]  made  a  significant  statement  in  his  article  about  polarimetric  properties 
in  general.  He  stated  that  if  the  target  and  the  clutter  had  differing  geometrical  characteristics, 
with  the  proper  polarimetric  processing,  the  radar  return  from  clutter  and  target  might  be  dis¬ 
tinguishable  [7 :pages  148-149].  One  algorithm  that  exploits  polarimetric  data  is  the  Polarimetric 
Whitening  Filter  (PWF).  More  in-depth  information  on  the  PWF  can  be  found  in  the  articles 
co-authored  by  Novak  [19,  20,  21].  These  processing  techniques  may  prove  useful  in  the  reduction 
of  speckle  within  the  SAR  images.  Novak  states  that  reduction  of  speckle  improves  the  ability 
to  detect  targets  in  the  image.  Novak  [21]  also  discusses  various  detection  algorithms  using  these 
various  processing  techniques.  He  discovered  that  these  processors,  the  Optimally  Weighted  Filter 
(OWF)  in  particular,  do  not  perform  significantly  better  than  working  with  a  single  polarization. 
He  did  state;  however,  that  a  “multi-look”  radar  system  should  improve  the  odds  of  detection. 
Since  WASAR  is  a  “multi-look”  system,  these  processing  techniques  could  prove  to  be  quite  useful. 

Target  detection  is  the  main  subject  of  this  thesis.  Detection  of  a  target  involves  the  deter¬ 
mination  of  the  presence  or  absence  of  the  target.  This  thesis  does  not  deal  with  the  identification 
of  a  target.  Since  we  don’t  know  any  specific  parameters  of  the  target,  the  problem  of  detection 
becomes  more  difficult.  However,  there  are  certain  assumptions  about  the  target  that  can  be  made. 
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By  assuming  a  target  is  man-made,  we  can  make  certain  speculations  about  its  characteristics. 
Man-made  objects  have  a  tendency  to  be  regularly  shaped,  that  is  composed  of  flat,  spherical, 
and/or  cylindrical  shapes.  Buildings  for  example  are  usually  composed  of  flat  walls  and  flat  or 
pyramidally  shaped  roofs.  Naturally  occurring  objects  usually  have  irregular  surfaces,  such  as  a 
forested  area,  a  shrub  filled  plain,  or  to  a  certain  extent  even  a  grassy  plain  is  relatively  irregular 
in  appearance  compared  to  a  building.  Smaller  man-made  objects,  such  as  vehicles,  also  appear  to 
be  regularly  shaped. 

Two  techniques  are  currently  used  for  target  detection;  change  detection  and  threshold  de¬ 
tection.  Change  detection  uses  two  separate  images  of  the  same  scene  to  determine  any  differences 
between  the  two  images.  Normally,  this  would  involve  a  temporal  change,  a  change  in  time  between 
the  two  images.  However,  it  is  possible  to  exploit  some  of  these  change  detection  techniques  when 
using  the  multiple  images  associated  with  WASAR.  Problems  arise  with  change  detection  when  the 
object  being  detected  doesn’t  vary  with  aspect  angle.  The  other  detection  technique,  threshold¬ 
ing,  involves  the  comparison  of  pixel  intensities  with  a  known  value  or  threshold.  All  pixels  above 
this  threshold  level  are  then  identified  as  target  pixels.  Problems  with  thresholding  techniques 
occur  when  the  intensities  of  the  clutter  and  target  are  indistinguishable  making  it  harder  to  set 
an  appropriate  threshold.  The  goal  of  this  research  is  to  use  a  combination  of  these  two 
techniques  in  tandem  with  an  optimal  polarimetric  processing  technique  to  perform 
target  detection. 

Difficulties  with  detection  in  general  using  WASAR  stem  from  the  relatively  recent  availability 
of  WASAR  data.  How  classical  detection  techniques  will  perform  on  WASAR  imagery  is  unknown. 
No  published  techniques  are  available  for  optimally  combining  the  WASAR  images.  Theoretically 
though,  clutter  returns  will  stay  relatively  the  same  throughout  aspect  changes;  whereas,  target 
returns  (e.g.  returns  from  flat  plates)  will  change  more  drastically  over  aspect  variations.  If 
this  information  could  be  exploited  by  using  some  type  of  change  detection  technique,  perhaps 
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target  detection  could  become  more  probable  than  using  a  single  image.  Though  most  techniques 
of  change  detection  [1,  6,  28]  involve  the  use  of  the  same  image  (same  frequency,  aspect  angle, 
depression  angle,  etc.  except  for  a  temporal  difference),  some  of  these  detection  techniques  may 
find  application  in  change  detection  between  aspect  angles. 

1.3  Concluding  Remarks 

As  stated  previously,  this  thesis  deals  with  exploiting  data  available  through  Wide-Angle  SAR 
for  the  purpose  of  target  detection.  The  initial  problem  is  data  availability,  since  Wide-Angle  SAR 
concepts  are  relatively  new.  In  order  to  work  with  WASAR  data,  the  data  must  first  be  available 
and  realistic.  The  second  part  of  the  problem  with  working  with  WASAR  images  reveals  itself 
logically  as  a  registration  problem  between  the  various  aspect  angles.  The  coordinate  system  in 
one  image  must  be  related  to  the  coordinate  system  in  another.  The  third  step  is  pre-processing 
in  which  the  images  are  combined  polarimetrically  to  optimize  any  following  detection  algorithm. 
This  step  is  an  attempt  to  use  the  polarimetric  diversity  advantageously  in  detecting  targets.  The 
last  step  of  course  is  the  detection  phase.  The  ultimate  goal  of  this  research  is  to  determine  if 
the  use  of  WASAR  instead  of  SAR  improves  the  ability  to  detect  targets.  These  four  areas  are 
discussed  in  the  following  chapters. 
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II.  Pre-Detection  Processing 


2.1  Data  Generation 

The  initial  step  in  the  detection  of  a  target  within  a  SAR  image  is  obtaining  a  SAR  image 
for  manipulation.  Since  real  WASAR  imagery  is  not  currently  available,  it  is  necessary  to  generate 
the  image  data  using  a  new  software  package  from  Loral.  The  software  allows  the  generation  of  an 
image  with  multiple  targets  (M35  trucks)  randomly  placed  throughout  the  image  region.  Various 
images  were  generated  for  detection  comparison  including:  targets  in  grass,  targets  in  a  scene 
with  25%  forestation,  and  targets  in  a  scene  with  50%  forestation.  Within  the  images,  any  forested 
region  contains  deciduous  trees  in  fall  season  with  variable  tree  density.  The  simulated  SAR  system 
operates  at  1200  MHz  (L-band)  with  constant  spatial  resolution.  The  output  is  fully  polarimetric 
complex  image  data.  An  L-band  frequency  is  used  because  this  range  of  frequencies  has  been  shown 
to  have  foliage  penetration  capabilities  over  the  use  of  higher  frequencies.  The  targets  are  randomly 
interspersed  throughout  the  entire  image.  Imagery  data  is  generated  for  a  depression  angle  of  15° 
at  aspect  angles  of  ±  45°,  ±  30°,  ±  15°,  and  0°.  This  data  is  then  input  to  a  MATLAB®  routine 
which  displays  the  log  magnitude  value  of  the  images  as  shown  in  Figure  2.1.  This  figure  shows  all 
21  images  involved  in  the  seven  aspect  fully  polarimetric  image  data  set.  More  specific  information 
on  the  WASAR  imagery  is  contained  in  Appendix  A. 

Within  the  scene,  each  pixel  represents  approximately  2  square  feet  of  the  actual  ground 
scene.  Each  image  is  composed  of  512  x  512  pixels.  The  +45°  HH  image  displayed  in  Figure  2.1 
is  also  shown  enlarged  in  Figure  2.9  on  page  2-15.  The  images  appear  to  show  a  forested  stretch 
of  land  within  a  grassy  area.  Notice  that  the  higher  returns  (whiter  areas)  correspond  to  the  trees 
and  targets. 
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2.2  Aspect  Angle  Registration,  Image  Rotation 


The  next  step  in  the  detection  process  is  to  devise  some  scheme  of  registering  the  various 
images  so  that  pixel  coordinates  in  one  aspect  relate  to  those  in  another.  If  the  images  are  not 
“ground  plane”  images,  this  could  be  a  very  difficult  problem  involving  projection  of  the  image 
into  the  ground  plane  and  interpolation.  Fortunately,  the  images  generated  for  this  thesis  are  in 
the  ground  plane  and  require  only  a  simple  rotation  in  order  to  register  pixel-to-pixel  between 
the  various  images.  This  assumption  was  verified  by  using  the  rotational  method  described  in  the 
following  paragraph  and  then  matching  objects  within  the  images.  It  may  also  be  of  interest  to 
note  that  the  center  of  rotation  necessary  to  match  the  images  turned  out  to  be  the  very  center  of 
each  image.  The  rotation  process  used  in  this  research  also  caused  the  image  registration  between 
various  aspect  angles  to  be  off  by  one  pixel  at  the  most.  This  means  that  an  object,  after  rotation 
of  the  image,  may  have  moved  to  a  position  that  was  shifted  by  one  pixel  from  its  true  position. 
This  offset  was  not  deemed  to  be  a  large  problem,  and  for  the  most  part  it  could  not  be  corrected 
since  each  image  is  a  discrete  set  of  pixels. 

The  images  generated  by  the  Loral  software  were  contained  in  a  512  x  512  matrix  of  pixel 
values.  The  choice  of  rotation  is  arbitrary,  as  long  as  all  aspects  are  rotated  to  the  same  coordinate 
system.  In  order  to  simplify  the  rotation  algorithm,  all  aspects  are  rotated  to  the  most  positive 
aspect,  45°.  This  means  that  one  image,  the  +45°  image,  does  not  have  to  be  rotated  at  all,  and 
the  —45°  image  can  be  rotated  using  ?o  intrinsic  MATLAB®  function.  After  rotation,  it  should  be 
obvious  to  the  reader  that  the  overlapping  or  common  area  between  the  different  image  aspects  is 
a  shape  approaching  a  circle.  Since  MATLAB®  wor*s  more  easily  with  square  images,  the  images 
are  reduced  to  a  square  area  that  fits  inside  of  this  overlapping  circular  region.  The  result  is  the 
previously  mentioned  512  x  512  images  are  now  reduced  to  363  x  363  pixel  images;  however,  the 
coordinates  for  one  image  map  directly  to  the  coordinates  in  any  other  image. 
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(a)  HH  +45®  (b)  HH  0®  (c)  HH  -45® 


(d)  HV  +45®  (e)  HV  0®  (f)  HV  -45® 


(g)  W  +45®  (h)  W  0®  (i)  VV  -45® 

Figure  2.2  Images  Displayed  After  Rotation 
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The  images  generated  from  the  last  section  are  shown  rotated  in  Figure  2.2.  Only  the  ±  45° 
and  0°  aspect  angles  are  shown.  Notice  that  the  images  are  very  similar  to  each  other;  however, 
there  are  unique  features  to  each  of  the  nine  images. 

2.3  Pre-processing 

After  the  images  are  rotated,  various  types  of  pre-processing  on  the  polarimetric  data  can  be 
done.  It  is  desirable  to  find  some  optimum  way  of  combining  the  HH,  HV,  and  VV  images  in  order  to 
optimize  target  detection.  Four  methods  of  combining  these  images  are  explored:  simple  average, 
span,  optimal  weighting,  and  polarimetric  whitening  filter.  The  span,  optimal  weighting,  and 
polarimetric  whitening  filter  were  all  proposed  by  Novak  [20].  For  derivations  of  these  particular 
pre-processing  techniques  please  refer  to  the  article  indicated.  The  Novak  techniques  are  used 
in  this  research  with  certain  qualifications:  1)  measured  statistics  were  used  instead  of  actual 
statistics,  and  2)  the  entire  image  was  used  and  was  not  separated  according  to  clutter  and  target 
when  determining  statistics.  Measured  statistics  were  used  instead  of  actual  data  mainly  because 
in  a  real-world  image,  all  that  is  available  is  the  measured  data.  Any  method  that  reduces  the 
background  clutter  and  enhances  the  target  return  is  beneficial  in  this  preprocessing  stage. 

2.3.1  Polarimetric  Average.  The  first  method  this  research  uses  is  the  polarimetric 
average,  which  is  simply  the  average  of  the  HH,  HV,  and  VV  images 

ave=|(|HH|  +  |HV|  +  |W|).  (2.1) 

Equation  2.1  represents  a  pixel-by-pixel  operation  carried  out  on  each  of  the  different  aspect  angles. 
The  |  •  |  symbol  indicates  a  pixel-by-pixel  magnitude,  or  modulus,  value  and  is  not  to  be  confused 
with  a  determinant  of  a  matrix.  The  images  that  result  from  this  method  Me  shown  in  Figure  2.3. 
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(a)  +45°  (b)  0°  (c)  -45° 

Figure  2.3  Images  Displayed  After  Rotation  and  Polarimetric  Averaging 

2.8.2  Polarimetric  Span.  The  next  method  explored  is  the  Polarimetric  Span,  which  is 
another  type  of  averaging  technique.  The  span  operation  is  given  by 

span  =  |HH|2  +  2  |HV|2  +  |W|2  .  (2.2) 

Due  to  the  squaring  within  Equation  2.2,  the  higher  intensity  values  are  emphasized  and  the  impact 
of  lower  intensity  returns  are  lessened.  Notice  that  the  results  from  a  span  operation  are  very  similar 
to  those  obtained  for  the  average  (see  Figure  2.4),  though  the  background  (grass)  is  darker,  implying 
that  the  contrast  has  been  increased.  Notice  that  in  liquation  2.2,  phase  information  is  completely 
lost.  If  the  actual  values  of  the  averaged  and  span-encoded  images  were  examined  and  compared, 
they  would  be  vastly  different  as  can  be  seen  by  Equations  2.1  and  2.2. 


(a)  +45° 


(c)  -45« 


Figure  2.4  Image  Displayed  After  Rotation  and  Polarimetric  Span 
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g.S.S  Polarimetric  Optimal  Weighting.  The  Polarimetric  Optimal  Weighting  function, 
takes  the  span  one  step  further  in  that  it  uses  statistical  information  from  the  image  itself  to  weight 
the  various  polarimetric  returns.  The  weighting  function  is  described  by 

pow  =  |HH|J  +  *2|HV|J  +  *3|VV|2.  (2-3) 


Where  the  constants,  and  is,  are  given  by 


l  +  |p|2 


*3=-, 

7 


E(HH.W') _  (26) 

{ |HH  |a  }  E  { |  VVI2} 

I-E(|HV|i}  (2.7) 

E{|HB|’}' 

-  E(|VV|’}  12.8) 

7  e{|HH|’}’ 

Notice  that  Equations  2.6,  2.7,  and  2.8  use  the  statistical  mean  of  the  image  (  E{-}  )-  This  is  a 
measured  mean  and  therefore  an  approximation  of  the  actual  mean  of  the  image.  As  the  data  set 
gets  larger,  approaches  infinity,  the  measured  mean  will  approach  the  actual  mean.  The  actual  mean 
would  be  the  mean  of  the  probability  density  function  which  was  used  to  create  the  cluttered  image. 
In  a  real-world  image  only  the  measured  data  is  available  for  determining  statistical  information; 
therefore,  it  is  more  prudent  to  test  algorithms  using  the  measured  statistics  even  if  the  actual 
statistics  are  available.  The  results  of  optimal  weighting  are  shown  in  Figure  2.5. 
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(a)  +45’  (b)  0°  (c)  -45° 

Figure  2.5  Image  Displayed  After  Rotation  and  Polarimetric  Optimal  Weighting 

2.3.4  Polarimetric  Whitening  Filter.  Finally,  the  Polarimetric  Whitening  Filter  is  de¬ 
scribed  by 


pwf  = 


\HH\* 


|VV|2 


\HV\ 2 


2|p| 


<thb(  1  -  M2)  0hh(  1  -  |p|2fr  0hh(  <thh(1  -  |p|2)>/7 


(2.9) 


where  p,  e,  and  7  are  as  defined  previously  in  Equations  2.6, 2.7,  and  2.8  respectively.  The  parameter 
<thh  i«  defined  as 

<thh  =  E{|HH|2}.  (2.10) 


The  parameters  <t>vv,  and  <t>p  are  the  phase  angles  of  the  indicated  subscripts.  The  results  of 
using  this  PWF  are  shown  in  Figure  2.6.  Notice  the  similarities  between  the  Optimal  Weighting 


(a)  +45°  (b)  0°  (c)  -45° 

Figure  2.6  Images  Displayed  After  Rotation  and  Polarimetric  Whitening 


2-8 


function  Equation  2.3  and  the  Polarimetric  Whitening  Filter  (PWF)  defined  by  Equation  2.9.  The 
PWF  takes  advantage  of  not  only  statistical  information  relating  to  the  image,  but  also  exploits 
phase  data.  The  images  as  they  are  displayed  in  this  figure  appear  to  be  very  similar  to  the 
optimally  weighted  images.  If  these  whitened  images  were  displayed  in  the  same  manner  as  the 
optimally  weighted  images,  the  display  would  be  almost  totally  white.  In  order  to  display  the 
whitened  images  and  be  able  to  see  any  detail,  the  intensity  of  the  display  was  diminished.  Most  of 
the  images  that  are  displayed  in  this  document  use  the  same  intensity  setting  so  that  comparisons 
can  be  made  between  the  various  images.  A  few  of  the  figures  do  use  a  different  image  intensity 
setting  or  imaging  method  though.  When  a  variation  in  the  imaging  technique  is  used,  a  note  will 
be  made  of  the  fact.  This  variation  in  intensity  does  not  affect  any  of  the  procedures  accomplished 
in  this  research  since  the  intensity  levels  are  adjusted  only  during  display  for  ease  of  viewing. 

2-4  Statistical  Analysis 

As  stated  previously,  the  polarimetric  combining  techniques  might  prove  useful  for  improving 
our  target  detection  capability.  In  order  to  measure  this  improvement,  some  metric  must  be  used  to 
compare  the  various  polarimetric  techniques.  To  this  end,  there  are  two  metrics  that  are  analyzed 
in  order  to  quantify  whether  the  previously  mentioned  polarimetric  techniques  improved  the  images 
in  some  way.  The  first  metric  is  the  standard  deviation  to  mean  ratio  of  the  image.  The  other 
metric  is  target  to  clutter  ratio.  The  target  to  clutter  ratio  is  important  since  any  decrease  in  this 
ratio  would  also  decrease  the  probability  of  detection. 

The  first  metric,  the  standard  deviation  to  mean  ratio,  is  defined  by  Novak  [20:page  294]. 
This  metric  is  computed  by  calculating  the  standard  deviation  of  the  pixel  intensities  of  the  entire 
image  and  dividing  by  the  mean  value  of  the  pixel  intensities  of  the  whole  image.  This  ratio  is  then 
converted  to  dB  and  used  to  compare  the  various  polarimetric  processing  techniques. 
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1ES23B 

+30° 

SI 

0° 

kb 

mmw 

Standard  Deviation  to  Mean  Ratio  (d 

B) 

HH 

41.42 

38.80 

35.28 

43.10 

46.34 

42.39 

39.15 

HV 

41.69 

35.69 

33.00 

38.05 

39.89 

33.27 

36.57 

VV 

43.75 

33.63 

32.77 

47.41 

39.95 

44.19 

45.57 

averaged 

36.18 

36.32 

33.54 

41.48 

37.17 

37.69 

37.54 

spaned 

30.46 

23.52 

23.43 

29.80 

30.75 

30.60 

30.44 

optimal 

28.14 

23.67 

23.73 

27.21 

28.24 

28.23 

28.05 

whitened 

27.08 

24.41 

24.38 

25.83 

27.05 

27.78 

27.02 

1 

target  to  Clutter  Ratio 

(dB) 

HH 

35.91 

38.12 

31.96 

40.43 

36.44 

34.27 

35.81 

HV 

32.93 

31.83 

31.57 

31.69 

33.56 

26.84 

29.44 

VV 

36.96 

30.76 

32.44 

32.07 

30.14 

36.41 

35.33 

averaged 

35.49 

35.51 

31.14 

35.88 

33.73 

35.35 

36.54 

spaned 

35.35 

32.42 

32.09 

35.81 

35.47 

34.88 

35.55 

optimal 

35.08 

33.05 

32.68 

35.31 

35.00 

34.66 

35.28 

whitened 

34.76 

33.26 

33.03 

34.89 

34.63 

34.55 

35.00 

Table  2.1  Statistics  for  Scene  1:  Targets  in  Grass 


The  second  metric,  target  to  clutter  ratio,  is  used  to  determine  if  any  of  the  polarimetric 
processing  techniques  enhance  target  detection  capabilities.  Since  the  target  locations  are  known, 
it  is  fairly  simple  to  find  the  intensity  of  the  pixels  within  a  square  of  11  x  11  pixels  centered  about 
the  known  target  coordinates.  An  11  x  11  square  was  chosen  to  insure  that  only  pixels  associated 
with  a  target  are  included  in  the  target  intensity  values.  The  largest  possible  square  was  chosen  in 
order  to  include  as  many  pixel  returns  as  possible  to  obtain  an  accurate  target  mean.  All  the  target 
pixel  intensities  are  averaged  for  a  particular  image  and  then  divided  by  the  average  pixel  intensity 
of  the  entire  image.  The  mean  of  the  image  is  taken  over  the  entire  image  (including  the  targets) 
with  the  assumption  that  the  image  is  luge  (in  number  of  pixels)  and  the  number  of  target  pixels 
is  small.  All  of  the  scenes  generated  for  this  research  contain  8  to  10  targets  and  the  image  size 
after  rotation  is  363  x  363.  Therefore,  since  the  number  of  targets  is  small  and  the  image  size  is 
relatively  large,  the  impact  of  the  target  returns  on  the  image  mean  is  minimal.  Also,  notice  from 
Table  2.1  that  the  target  to  clutter  ratios  are  very  large  anyway.  Subtraction  of  the  targets  from 
the  clutter  mean  would  only  serve  to  increase  the  ratio  even  higher. 
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Figure  2.7  Scene  1:  Targets  in  Grass  (+45°  Aspect  HH  Image  After  Rotation) 


2-4.1  Scene  1:  Targets  in  Grass.  The  first  scene  analyzed  contains  targets  in  grass.  The 
+45°  HH  image  is  shown  in  Figure  2.7.  The  target  locations  are  designated  in  the  figure  using 
Xs.  The  returns  from  the  trucks  are  slightly  to  the  left  and  center  of  each  X.  Coordinates  of  the 
targets  are  available  in  a  file  generated  by  the  Loral  software.  Notice  that  there  appears  to  be  a 
return  without  an  X  beside  it.  This  is  most  likely  a  stray  tree  that  was  placed  in  the  scene  by  the 
generation  software.  The  software  designer  was  contacted  about  this  problem  and  the  existence  of 
a  tree  within  an  all  grass  scene  could  not  be  confirmed;  however,  the  designer  did  state  that  it  was 
possible. 

The  standard  deviation  to  mean  and  target  to  clutter  ratios  for  scene  1  are  found  in  Table  2.1. 
Notice  that  in  the  standard  deviation  to  mean  ratio  portion  of  the  table,  the  relative  values  in 
the  table  decrease  as  you  go  down  the  columns  of  the  processed  images.  This  result  is  expected 
since  the  polarimetric  processing  techniques  are  designed  to  reduce  this  metric.  As  expected,  the 
pclarimetric  whitening  filter  performs  the  best  on  average  in  reducing  standard  deviation  to  mean 
ratios.  Additionally,  target  to  clutter  ratio,  stays  relatively  the  same  (or  slightly  decreases)  down 
the  columns  of  the  Table  2.1.  This  trend  indicates  that  the  polarimetric  processing  techniques  do 
not  improve  target  detection  chances  appreciably,  and  in  some  cases  may  actually  decrease  the 
odds.  A  decrease  in  the  target  to  clutter  ratio,  results  in  the  target  becoming  more  and  more 
indistinguishable  from  clutter  intensities. 

2.4-2  Scene  2:  Targets  in  Trees  (25%  foliage).  The  second  scene  analyzed  contains  targets 
in  trees  and  grass  with  25%  forestation.  The  +45°  HH  image  is  shown  in  Figure  2.8.  Table  2.2 
indicates  that  the  techniques  proposed  by  Novak  [20]  -  span,  optimal  weighting,  and  pwf  -  perform 
the  best  on  this  data  set  for  reducing  standard  deviation  to  mean  ratios.  The  target  to  clutter  ratio 
decreases  slightly  from  the  un-processed  image,  and  in  this  case,  the  averaging  technique  performed 
the  best. 
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Figure  2.8  Scene  2:  Targets  in  Trees  and  Grass  (25%  foliage)  (4-45°  Aspect  HH  Image  After 
Rotation) 
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HH 

38.24 

40.31 

46.33 

39.68 

38.33 

30.80 

34.84 

HV 

33.88 

30.79 

33.67 

27.10 

27.82 

29.86 

32.56 

VV 

38.60 

34.25 

31.49 

36.93 

45.70 

30.68 

32.17 

averaged 

35.20 

45.24 

39.56 

35.04 

34.26 

28.11 

30.06 

spaned 

17.47 

14.63 

19.43 

26.52 

26.60 

18.36 

19.64 

optimal 

17.55 

17.01 

18.17 

23.93 

24.50 

19.05 

19.12 

whitened 

17.89 

17.51 

18.69 

22.69 

24.25 

18.45 

19.59 

1 

?arget  to  Clutter  Ratio 

(dB) 

HH 

37.74 

37.27 

46.92 

39.01 

35.82 

33.06 

34.27 

HV 

16.07 

30.20 

33.32 

28.83 

24.84 

30.05 

25.43 

VV 

40.14 

37.35 

32.78 

20.08 

45.75 

32.03 

29.41 

averaged 

39.38 

46.26 

37.48 

34.42 

37.95 

31.69 

31.77 

spaned 

30.05 

27.55 

28.64 

35.71 

33.50 

29.78 

30.87 

optimal 

30.74 

29.55 

28.88 

34.31 

33.24 

30.90 

31.71 

whitened 

30.89 

29.61 

29.10 

33.84 

33.09 

30.75 

31.84 

Table  2.2  Statistics  for  Scene  2:  Targets  in  25%  Foliage 


Notice  that  on  average,  the  target  to  clutter  ratio  has  dropped  as  compared  to  targets  in 
grass  alone  (compare  Tables  2.1  and  2.2).  This  result  is  expected.  It  is  reasonable  to  believe  that 
a  target  in  the  forest  would  be  a  little  bit  harder  to  detect  than  that  same  target  in  grass.  Also 
notice  that  the  standard  deviation  to  mean  ratios  are  lower  for  scene  1  than  for  scene  2. 


2-4-3  Scene  3:  Targets  in  Trees  (50%  foliage).  The  third  scene  analyzed  contains  targets 
in  trees  and  grass  with  50%  forestation.  The  +45®  HH  image  is  shown  in  figure  2.9.  The  statistical 
information  is  contained  in  Table  2.3.  Information  is  missing  from  the  +30®  aspect  due  to  the 
corruption  of  the  VV  polarized  image  (see  the  +30®  W  image  displayed  in  Figure  2.1  on  page  2- 
2).  The  cause  of  this  corruption  is  unknown,  but,  was  a  function  of  the  parameters  used  in  the 
generation  of  this  particular  scene.  It  may  also  be  of  interest  to  the  reader,  this  is  the  scene  that 
is  used  for  the  figures  contained  in  the  pre-processing  section. 

Notice  in  Table  2.3  that  the  standard  deviation  to  mean  ratio  experiences  large  reductions  as 
the  pre-processing  technique  becomes  more  complex,  moving  down  the  columns.  The  polarimetric 
whitening  filter  again  offers  the  best  performance  for  this  metric.  Also,  notice  that  the  target  to 
clutter  ratio  is  not  reduced  significantly  by  the  pre-processing  stage.  However,  comparing  Tables  2.1, 
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Figure  2.9  Scene  3:  Target  in  Trees  and  Grass  (50%  foliage)  (+45°  Aspect  HH  Image  After 
Rotation) 
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spaned 


optimal 


whitened 


HH 


HV 


VV 


averaged 


spaned 


optimal 


whitened 


Standard  Deviation  to  Mean  Ratio 


24.37 


28.43 


30.46 


22.82 


10.19 


12.88 


13.23 


Target  to  Clutter  Ratio 


24.94  -  25.24  41.28 


24.36  -  25.00 


31.43  -  32.12 


25.48  -  18.97 


22.84  -  23.22 


24.43  -  23.33 


24.80  -  23.79 


25.48 

39.41 

33.57 

25.66 

33.45 

26.27 

27.96 

24.77 

27.32 

23.19 

35.03 

33.87 

31.71 

23.90 

32.25 

22.82 

28.47 

30.09 

21.84 

25.01 

17.03 

25.13 

13.26 

12.03 

11.65 

15.80 

23.07 

15.24 

13.25 

13.24 

16.66 

22.14 

14.96 

12.53 

13.84 

31.95 

26.30 

37.47 

28.06 

26.92 

22.06 

37.44 

25.85 

21.40 

36.27 

22.46 

25.49 

24.51 

22.98 

24.30 

26.17 

24.55 

25.75 

26.17 

24.57 

26.10 

Table  2.3  Statistics  for  Scene  3:  Targets  in  50%  Foliage 


2.2,  and  2.3,  we  can  see  a  large  drop  in  target  to  clutter  ratios.  This  is  to  be  expected  since  the 
forested  portion  of  our  scene  is  becoming  thicker.  If  the  entire  scene  were  covered  with  forest,  the 
target  to  clutter  ratio  might  be  expected  to  drop  even  further,  perhaps  down  to  unity.  If  this  were 
to  happen,  say  targets  in  a  very  thick  forest,  detection  would  become  very  difficult.  This  extreme 
case  could  not  be  explored  since  the  Loral  software  only  allowed  a  maximum  of  50%  forestation. 

2.4  4  Statistical  Overview.  The  purpose  of  any  polarimetric  pre-processing  is  to  improve 
target  detection  capability  (over  the  use  of  a  single  polarimetric  image).  It  would  appear  from 
the  three  data  sets  that  were  analyzed  in  this  research  that  the  polarimetric  pre-processing  for  the 
most  part  does  not  appreciably  affect  target  detection  capability  as  quantified  by  the  target  to 
clutter  ratios.  The  use  of  a  pre- preceding  stage  does  not  improve  the  target  to  clutter  ratio,  and 
in  some  cases  can  actually  decrease  the  ratio.  However,  the  decreases  resulting  from  this  research 
are  not  unacceptable.  A  decrease  would  be  considered  unacceptable  if  the  target  to  clutter  ratio 
were  to  approach  unity.  The  results  from  this  chapter  are  inconclusive  in  favoring  one  technique 
over  another  for  optimization  of  target  detection  since  the  target  to  clutter  ratio  stays  relatively 
constant. 
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III.  Target  Detection 


The  main  purpose  of  this  thesis  is  target  detection.  To  this  end,  threshold  detection  is 
accomplished  on  individual  aspect  angles.  Threshold  detection  in  its  basic  form  is  simply  the 
comparison  of  the  pixel  intensity  levels  with  a  preset  value  or  threshold.  If  the  intensity  of  the 
pixel  being  tested  is  above  this  threshold,  then  that  pixel  is  returned  as  a  “target”.  After  the 
threshold  detector  is  used  on  all  the  image  aspects,  a  comparison  of  the  returned  pixel  coordinates 
is  accomplished.  For  this  thesis,  this  comparison  was  simply  an  accumulation  of  all  the  coordinates. 
In  set  terminology,  this  is  known  as  a  union  set  of  the  coordinates  of  the  individual  aspects.  The 
threshold,  7,  is  set  by  arbitrarily  selecting  a  value  defined  between  the  maximum  and  minimum 
pixel  intensities 

7  =  min {|pixel  value|}  +  [maxflpixel  value|}  -  min{|pixel  value|}]a.  (3.1) 

where  a  takes  on  values  ranging  from  0  to  1.  In  the  results  to  follow,  the  parameter  a  is  plotted  as 
the  independent  variable.  The  parameter  a  is  referred  to  as  the  threshold  and  allowed  to  range  from 
zero  to  unity.  Notice  that  this  equation  and  the  realization  in  the  code,  “detect2.m”  on  page  B-47, 
is  different.  The  minimum  pixel  values  range  from  10-5  to  10~3;  whereas,  the  maximum  is  generally 
in  the  102  range.  Which  indicates  that  the  minimum  is  approximately  zero  when  compared  to  the 
maximum  pixel  value  and  can  be  ignored  in  the  formula,  making  the  code  realization  equivalent  to 
Equation  3.1. 

Two  detection  methods  are  used.  For  the  first  method  of  detection,  the  threshold  is  found 
for  the  single  +45°  aspect  image  and  then  used  for  detection  on  the  other  aspects.  In  essence,  this 
procedure  sets  a  permanent  threshold  at  some  absolute  value.  For  the  second  method  of  detection, 
a  separate  threshold  is  derived  using  Equation  3.1  for  each  image.  This  second  method,  then, 
allows  a  “floating”  threshold  of  sorts.  Once  again,  detection  is  done  on  individual  aspects  and  the 
target  coordinates  combined.  Threshold  detection  is  also  accomplished  on  one  image  aspect  (the 
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+45°  aspect  in  particular)  and  the  results  compared  to  detection  using  all  the  image  aspects,  which 
equates  to  a  comparison  of  SAE  versus  WASAR. 

Also,  it  may  be  of  interest  to  note,  a  target  is  considered  detected  if  the  detection  algorithm 
returns  a  coordinate  within  IS  pixels  of  the  known  coordinate.  Since  the  coordinates  of  the  targets 
are  known,  it  is  possible  to  determine  the  pixel  size  of  the  returns  associated  with  the  targets. 
Several  of  the  targets  were  observed  in  several  aspect  angles  and  they,  the  targets,  were  found  to 
have  returns  that  ranged  in  size  from  12x12  up  to  15x15  pixels.  So  a  pixel  patch  of  20x20  was 
chosen  to  ensure  that  any  pixels  within  this  block  would  be  grouped  into  the  same  target.  Of 
course,  this  pixel  patch  would  need  to  change  with  different  size  targets.  The  size  of  this  patch  also 
limits  the  target  resolution  capability;  the  ability  to  distinguish  between  two  targets  that  are  close 
together. 

The  four  pre-processing  techniques  discussed  in  the  last  chapter  are  used  as  well  as  simply 
using  the  HH  polarized  image  by  itself.  Detection  after  using  these  various  pre-processing  stages  also 
determines  whether  any  of  the  pre-processing  techniques  improve  the  target  detection  capability. 
Any  large  variations  in  the  detection  performance  should  stand  out  if  the  threshold  value  is  varied 
over  a  wide  range. 

Figure  3.1  shows  an  enlarged  view  of  100  x  100  pixel  area  around  a  typical  target  within  the 
images  generated  in  this  research.  The  target  is  the  patch  of  high  intensity  returns  in  the  center 
of  the  image.  As  stated  previously  this  is  an  M-35  truck  and  there  can  be  as  many  as  15  trucks  in 
a  scene  at  various  target  aspect  angles  to  the  probing  radar  system.  Notice  that  there  are  other 
pixels  within  the  image  that  are  very  close  in  intensity  to  the  target  pixels.  Herein  lies  the  crux  of 
the  target  detection  problem.  If  clutter  pixels  have  intensities  as  high  as  the  target’s,  then  a  simple 
threshold  detector  identifies  these  clutter  pixels  as  targets,  resulting  in  false  alarms. 
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Figure  3.1  Enlargement  of  Typical  Target 
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S.l  Scene  1:  Targets  in  Grass 


Figures  3.2  and  3.3  show  the  results  of  using  multiple  aspect  angle  detection.  The  figures  show 
five  different  plots  which  indicate  the  results  of  using  various  pre-processing  techniques.  Notice  that 
for  the  first  scene,  8  targets  are  common  to  all  of  the  rotated  image  sets.  Also,  notice  that  when 
targets  are  out  in  the  open  as  they  are  in  this  scene,  false  alarms  are  not  a  big  problem.  When  man¬ 
made  objects  are  placed  in  grass  (or  shrubs  for  that  matter),  setting  the  threshold  is  not  a  major 
problem.  Notice  that  false  alarms  do  not  occur  in  this  algorithm  until  the  threshold  is  set  below 
0.2.  At  this  point  the  threshold  is  becoming  lower  than  the  intensity  values  present  in  the  grassy 
area.  With  this  obviously  easy  detection  situation,  the  SAR  and  WASAR  methods  of  detection 
work  relatively  the  same,  except  that  with  WASAR,  the  threshold  can  operate  over  a  much  wider 
range  and  still  detect  all  of  the  targets.  Notice  that  both  the  SAR  and  WASAR  are  able  at  some 
point  to  detect  all  the  targets  and  at  the  same  time  not  realize  any  false  alarms.  This  point  in  the 
detection  and  false  alarm  curves  will  not  always  be  an  easy  situation  to  attain,  as  will  be  observed 
in  scenes  two  and  three  to  follow.  The  two  detection  techniques,  pre-set  in  Figure  3.2  and  scaled 
threshold  in  Figure  3.3,  are  fairly  indistinguishable  in  detection  performance.  Both  techniques  are 
able  to  attain  a  point  at  which  all  the  known  targets  are  detected  and  no  false  alarms  are  indicated. 

A  particular  anomaly  of  interest  can  be  found  in  the  Averaged  plot  in  Figures  3.2.  Notice  that 
between  0.5  and  0.6  threshold  values,  the  number  of  targets  detected  is  7;  whereas,  at  points  greater 
and  less  than  these  thresholds  the  number  of  targets  detected  is  8.  This  would  seem  intuitively 
impossible.  It  indicates  that  as  the  threshold  is  lowered  at  this  point,  instead  of  increasing  the  odds 
of  detecting  a  target,  a  target  is  lost.  This  can  be  easily  explained.  As  pixels  are  grouped  together 
to  form  a  target  after  pixel  detection,  their  coordinates  are  simply  averaged  together.  The  grouping 
is  within  a  20  pixel  block  as  explained  previously.  Now,  if  two  targets  happen  to  be  very  close, 
say  within  the  20  pixel  block  of  each  other  or  touching,  these  two  targets  will  get  grouped  into 
one  coordinate.  This  is  what  is  known  as  target  resolution,  the  ability  to  distinguish  between  two 
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Figure  3.2  Detection  Results  for  Targets  in  grass  (pre-set  threshold) 
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Figure  3.3  Detection  Results  for  Targets  in  grass  (scaled  threshold) 
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closely  spaced  targets.  This  definitely  is  a  problem  for  a  real  world  system,  since,  ideally,  all  targets 
should  be  identified  and  associated  coordinates  placed  near  the  center  of  the  target.  However,  for 
the  purposes  of  this  thesis,  investigating  the  utility  of  WASAR  over  SAR,  this  problem  is  not  fated. 
All  we  need  do  is  take  note  of  the  anomaly  and  ensure  that  the  cause  does  not  affect  the  algorithm, 
which  it  doesn’t. 

3.2  Scene  2:  Targets  in  Trees  and  Grass  (25%  foliage) 

Figure  3.4  shows  an  enlarged  view  of  100  x  100  pixel  area  around  a  typical  tree  within  the 
images  generated  in  this  research.  The  tree  is  the  patch  of  high  intensity  returns  in  the  center  of 
the  image.  All  of  the  trees  in  the  images  in  this  research  are  deciduous  trees  in  fall  season  with  a 
varying  trunk  density  between  trees.  The  number  of  trees  in  the  scene  depends  on  the  percentage 
of  foliation  selected  during  the  generation  of  the  image.  Notice  the  similarity  between  the  return  for 
a  tree,  Figure  3.4,  and  that  of  a  target,  Figure  3.1.  Both  seem  to  be  of  the  same  intensity,  though 
of  different  shapes.  This  can  be  a  major  problem  in  detecting  targets  buried  in  forest  clutter.  If 
the  target  does  not  have  a  considerably  higher  return  than  the  trees,  it  can  be  very  difficult  to  run 
a  simple  threshold  detection  technique  over  the  image.  On  the  fringe  areas  of  Figure  3.4  more  trees 
can  be  seen  in  the  area.  It  may  be  of  particular  interest  to  note,  the  high  intensity  return  on  the 
upper  edge  near  the  left  is  a  tree  and  a  truck  close  together. 

Figures  3.5  and  3.6  shows  the  results  of  using  multiple  aspect  angle  detection.  Nine  targets 
are  known  to  be  ^resent  in  this  scene.  Threshold  values  are  determined  for  the  cross-over  point, 
defined  as  the  point  at  which  the  number  of  targets  detected  equals  the  number  of  false  alarms.  This 
point  is  important  since  at  the  cross-over,  anything  detected  has  a  50/50  chance  of  actually  being 
a  target.  The  cross-over  point  can  be  clearly  seen  in  all  of  the  plots  for  both  the  SAR  and  WASAR 
processed  images.  Notice  that  the  false  alarm  rate,  the  rate  at  which  the  number  of  false  alarms 
increases  as  the  threshold  value  is  decreased,  seems  to  be  the  same  for  both  SAR  and  WASAR 
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Figure  3.5  Detection  Results  for  Targets  in  25%  foliage  (pre-set  threshold) 
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Figure  3.6  Detection  Results  for  Targets  in  25%  foliage  (scaled  threshold) 
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even  though  for  WASAR  the  false  alarms  start  at  a  higher  threshold.  The  most  important  thing 
to  notice  is  that  at  the  higher  thresholds,  past  the  point  at  which  false  alarms  occur,  the  WASAR 
is  able  to  detect  more  of  the  targets  on  average  than  the  SAR.  This  would  support  the  hypothesis 
that  the  threshold  can  be  set  higher,  to  decrease  false  alarms,  and  still  pick  up  most  of  the  targets. 
These  plots  show  that  at  some  point  within  the  various  image  aspects,  most  of  the  targets  will 
“flare-up”  higher  than  the  surrounding  clutter.  So,  with  WASAR  a  few  targets  can  be  missed  in  a 
single  aspect  image  and  then  picked  up  in  another  aspect. 

Notice  the  difference  in  performance  of  the  two  detection  methods.  The  first  item  to  notice 
is  that  with  scaled  thresholding,  the  number  of  targets  detected  goes  to  zero  when  the  threshold 
is  raised  to  1  (100%).  This  is  to  be  expected  on  a  single  aspect  thresholding  system  such  as  the 
scaled  thresholding  technique.  As  the  threshold  is  raised  toward  the  maximum  pixel  value,  the 
number  of  targets  detected  drops  off.  The  significant  comparison  on  this  point  however  is  observed 
by  using  the  preset  thresholding  technique  charts.  As  can  be  clearly  seen  in  Figure  3.5,  the  number 
of  targets  detected  using  WASAR  does  not  drop  to  zero  at  100%  cut-off.  This  would  indicate  that 
in  at  least  one  of  the  images,  the  target  pixel  intensity  (in  absolute  measurement)  is  higher  than  the 
maximum  pixel  intensity  in  the  +45°  image.  Remember  that  in  the  preset  thresholding  technique 
the  threshold  is  set  using  only  the  most  positive  aspect  image  in  the  image  set.  Also,  notice  that 
the  various  processing  techniques  yield  varying  results.  The  optimal  weighting  and  the  whitening 
filter  both  pull  the  false  alarm  curve  zero  cross-over  to  the  right.  For  the  start  of  the  false  alarm 
curves,  the  Spanned  images  seem  to  perform  the  best.  The  point  at  which  the  first  false  alarm 
occurs  is  quite  low,  at  a  threshold  of  0.2. 

Notice  the  differences  in  the  performance  of  the  detection  algorithm  on  scene  one,  Figures  3.2 
and  3.3,  as  compared  to  the  performance  on  scene  two,  Figures  3.5  and  3.6.  As  the  number  of  trees 
in  the  scene  increase,  the  detection  of  targets  becomes  more  difficult.  This  is  not  an  unexpected 
result.  As  the  targets  become  more  and  more  obscure,  it  will  be  more  difficult  to  pick  them  out. 
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S.S  Scent  S:  Targets  in  Trees  and  Grass  (50%  foliage) 


Figures  3.7  and  3-8  show  the  results  of  using  multiple  aspect  angle  detection.  10  targets  are 
common  to  all  of  the  image  sets.  Notice  the  difficulty  with  choosing  an  optimal  threshold  for  this 
scene.  At  no  point  can  we  attain  the  ideal  situation  in  which  all  the  known  targets  are  detected 
and  no  false  alarms  are  indicated.  Notice  in  Figure  3.7  that  the  Spanned  case  comes  close.  If  the 
cutoff  is  set  at  30%  for  a  scaled  threshold,  all  ten  targets  are  detected  and  one  false  alarm  results. 
A  few  of  the  other  cases  come  close,  but  none  realize  the  important  goal  of  no  false  alarms  and 
all  targets  detected.  It  is  also  interesting  to  note,  that  the  WASAR  seems  to  perform  better  than 
using  normal  SAR.  WASAR  comes  closer  in  all  cases  to  detecting  all  of  the  targets  and  except  for 
a  few  cases,  the  false  alarm  rates  appear  to  be  similar;  though,  once  again,  WASAR  false  alarms 
begin  at  a  higher  threshold. 

Notice  the  similarity  of  the  detection  performance  between  scenes  two  and  three,  Figures  3.5, 
3.6,  3.7,  and  3.8.  Both  are  very  similar  in  the  rates  at  which  detection  falls  off  and  the  rate  at 
which  false  alarms  increase.  This  similarity  would  indicate  that  the  amount  of  forestation  does  not 
impact  the  detection  performance  of  these  two  techniques  appreciably.  Specific  points  of  interest 
(points  at  which  false  alarms  drop  to  zero,  the  cross-over  point,  and  point  at  which  all  targets  are 
detected)  do  differ  though.  As  expected,  when  more  trees  are  present,  the  false  alarms  start  at 
higher  cut-offs,  the  point  at  which  all  targets  are  detected  starts  at  a  lower  cut-off,  and  so  on. 

The  results  from  the  previous  sections  indicate  that  WASAR  definitely  gives  better  perfor¬ 
mance  over  the  use  of  single-aspect  SAR.  The  actual  performance  of  the  detection  algorithm  is 
dependent  on  the  amount  of  forestation  present  in  an  image.  Also,  the  performance  of  the  pre-set 
and  scaled  methods  are  fairly  indistinguishable. 
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Figure  3.7  Detection  Results  for  Targets  in  50%  foliage  (pre-set  threshold) 
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Figure  3.8  Detection  Results  for  Targets  in  50%  foliage  (scaled  threshold) 
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3.4  Subtraction  Detection  Method 

Another  detection  technique  is  explored  but  the  results  are  inconclusive.  This  other  technique 
involves  the  comparison  of  the  various  aspect  angles  within  WASAE  by  subtraction.  In  this  research, 
the  method  of  subtraction  explored  is  one  in  which  all  aspects  are  subtracted  individually  from  the 
most  positive  aspect  image.  The  theory  behind  subtraction  is  that  across  the  aspects,  the  returns 
of  the  clutter  will  not  vary  appreciably;  whereas,  the  target  returns  will  vary  considerably  from 
one  aspect  to  another.  A  simple  understanding  of  this  approach  can  be  gained  from  an  example 
using  the  visible  light  spectrum.  No  matter  from  which  aspect  angle  a  tree  is  observed,  the  tree 
looks  very  similar.  Now  this  is  not  to  say  that  there  aren’t  variations  in  the  angular  views.  The 
cross-sectional  size  of  the  tree  may  change,  the  apparent  extent  of  the  foliage  coverage  may  change, 
etc.  However,  generally  speaking,  the  tree  looks  relatively  the  same.  Now  compare  this  to  various 
angular  views  of  say  a  semi-truck.  The  trucks  cross-sectional  size  changes  greatly  from  say  a  head 
on  view  compared  to  a  side  view.  If  the  truck  is  a  flat  bed,  it  may  appear  on  the  average  to  be 
very  tall  from  a  front  view  but  much  lower  on  average  from  the  side. 

Now  if  we  extend  this  example  into  the  realm  of  radar  waves,  we  can  theorize  as  to  the 
variation  of  the  tree  and  truck  while  being  illuminated  by  a  radar  beam.  The  main  return  from 
a  tree  is  the  result  of  the  radar  wave  bouncing  off  of  the  tree  trunk.  In  a  simple  sense,  a  tree 
trunk  is  a  cylindrical  object.  If  we  illuminate  this  cylindrical  object  from  various  aspect  angles 
perpendicular  to  the  cylindrical  axis,  the  return  looks  the  same.  The  next  highest  return  on  trees 
are  the  result  of  the  radar  wave  bouncing  off  of  various  sized  and  oriented  branches.  If  we  assume  a 
random  pattern  to  these  branches,  the  conglomerate  return  of  the  combined  branches  will  appear 
to  be  somewhat  of  a  spherical  shape,  or  conical  in  the  case  of  coniferous  trees.  Either  one  of  these 
shapes  will  clearly  yield  similar  returns  across  various  aspect  angles.  The  last  variable  to  account 
for  is  the  leaves.  The  leaves  will  have  a  tendency  to  yield  a  much  smaller  return,  but  still  be  similar 
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to  the  branches  in  conglomerate  shape  and  size.  Therefore,  the  return  here  is  also  similar  across 
the  aspect  angles.  So,  a  tree  theoretically  will  have  very  little  variation  across  aspect  angles. 

The  truck  on  the  other  hand,  when  illuminated  by  the  radar  will  vary  greatly  with  aspect 
angle.  We  can  assume  that  the  truck  is  composed  mainly  of  flat  plates  and  we  know  from  radar 
measurements  that  a  flat  plate  yields  its  highest  return  for  radar  illuminations  that  are  perpendic¬ 
ular  to  the  surface.  So,  the  truck  will  have  high  returns  at  the  sides  and  fronts  and  smaller  returns 
result  when  illuminated  at  skew  angles  of  the  flat  sides  and  front  -  when  the  radar  is  coming  in  at 
a  comer  of  the  truck.  Therefore,  the  truck  has  largely  varying  aspect  returns. 

Now  if  we  have  a  scene  in  which  natural  clutter  and  trucks  are  present  and  we  view  it  from 
various  angles,  the  clutter  will  appear  to  be  similar  in  all  scenes.  However,  the  truck  return  will 
change  from  angle  to  angle.  Theoretically  if  these  two  scenes  were  to  be  subtracted,  the  clutter 
return  should  go  to  zero  and  the  target  should  remain  relatively  high.  In  the  situation  in  this 
research  in  which  the  coordinates  of  the  two  images  may  be  off  by  a  pixel  due  to  the  rotation 
method  used,  the  subtraction  may  not  be  as  effective.  However,  if  there  is  a  mismatch  of  pixels, 
the  bulk  of  each  tree  should  still  subtract  with  the  possibility  of  a  high  intensity  ring  around  the 
tree  position  that  was  not  subtracted  away. 

Statistically  speaking,  subtraction  is  similar  to  a  summing  operation.  In  other  words,  the 
new  mean  of  the  subtracted  images  will  be  the  difference  of  the  individual  means.  If  we  assume 
that  the  two  clutter  scenes  are  identically  distributed  random  variables  then  the  mean  becomes 
zero.  The  variance  of  the  sum  of  two  identically  distributed  random  variables  is  twice  the  single 
image  variance.  So,  this  tells  us  that  the  clutter  variation  has  become  larger.  Now  the  targets  when 
subtracted  will  yield  the  difference  between  the  two  aspect  angles.  Difficulties  at  this  point  may 
arise  if  the  difference  in  the  target  aspects  is  not  greater  than  the  variation  in  the  clutter.  In  this 
case  the  subtracted  target  is  being  swamped  by  the  new  clutter  variation. 
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This  subtraction  technique  is  used  for  this  research  with  only  a  slight  variation.  After  sub¬ 
traction  of  an  individual  image  from  the  most  positive  image,  a  minimum  pixel  value  is  added  back 
into  the  image  to  ensure  that  no  zero  value  pixels  are  present.  This  is  done  purely  for  imaging 
purposes  and  in  no  way  impacts  on  the  usability  of  the  results.  After  subtraction  and  addition 
of  the  minimum,  the  resultant  image  is  run  through  a  threshold  detector.  The  detection  method 
described  earlier  in  this  chapter  as  the  scaled  threshold  method  is  used.  This  subtraction  tech¬ 
nique  is  explored  only  to  see  if  any  advantages  existed  over  the  more  simple  cumulative  detection 
techniques  discussed  in  the  previous  sections. 

For  this  research,  the  subtraction  method  is  used  on  scene  two  since  it  is  somewhat  of  an 
average  scene  of  the  three  generated.  The  results  of  subtraction  on  the  HH  rotated  image  data 
set  are  shown  in  Figure  3.9.  The  image  in  the  center  at  the  top  is  the  original  +45°  HH  image. 
All  the  images  in  this  figure  are  displayed  using  an  auto-scaled  imaging  technique.  This  means 
that  the  image  intensities  are  automatically  scaled  to  the  highest  and  lowest  pixel  intensities  to 
insure  that  all  of  the  pixel  values  are  displayed.  This  was  necessary  since  after  subtraction,  the 
pixels  will  have  a  tendency  to  be  grouped  very  tightly  around  a  few  intensity  values  with  very  little 
variation,  as  compared  to  an  un-subtr acted  image.  If  the  images  were  displayed  using  a  simple 
intensity  adjustment,  the  images  would  show  very  little  detail  since  they  would  appear  to  be  a 
uniform  intensity  across  the  entire  image. 

Notice  in  Figure  3.9,  that  the  contrast  in  the  images  is  very  small.  This  small  contrast 
indicates  that  our  hypothesis  on  subtraction  is  correct,  the  clutter  is  subtracting  from  itself.  In  a 
few  of  the  subtractions,  the  targets  really  stand  out.  Notice  especially  in  the  +45°— (0°)  subtraction, 
three  of  the  targets  can  be  clearly  seen  as  high  intensity  pixels  in  three  of  the  corners.  In  order  to 
pick  out  the  target  locations  in  the  darker  images,  compare  the  lighter  images  where  the  “Xs”  can 
be  clearly  seen  to  the  darker  images.  Notice  also,  that  in  some  of  the  subtracted  images,  the  trees 
do  still  have  return  values,  +45°  —  (+15°)  and  +45°  -  (-45°). 
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HH:  +45°  -  (+30°) 


HH:  +45®  -  (+15°) 


HH:  +45®  -  (0®) 


Figure  3.9  Images  Displayed  After  subtraction 
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(c)  Spanned 


TMlSH 

(d)  Optimal 


=====  9  Known  Targets 

- WASAR  Detected 

. WASAR  False  Alarms 

- SAR  Detected 

. SAR  False  Alarms 


(e)  Whitened 


Figure  3.10  Detection  Results  for  Targets  in  25%  foliage  (subtraction) 
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Figure  3.10  indicates  the  results  of  the  subtraction  detection  method.  If  these  results  are 
compared  to  those  found  in  Figure  3.6  you  will  notice  that  the  threshold  values  for  subtraction  are 
higher;  meaning  that  for  WASAR,  for  the  same  number  of  targets  detected,  the  threshold  must 
be  set  higher.  This  higher  threshold  value  also  occurs  for  the  number  of  false  alarms.  Notice  that 
the  SAR  curves,  both  number  detected  and  number  of  false  alarms,  are  exactly  the  same  in  both 
Figure  3.6  and  3.10.  This  is  an  expected  result  since  for  SAR  only  one  image  is  available  and 
no  others  are  available  for  subtraction.  Therefore,  for  the  single-aspect  SAR  case,  the  method  of 
detection  is  exactly  the  same  as  that  used  in  the  previous  sections,  and  the  results  are  exactly  the 
same.  The  results  in  Figure  3.10  indicate  that  for  this  particular  method  of  threshold  detection, 
subtraction  does  not  appreciably  improve  our  odds  of  detection.  However,  these  results  do  not 
exclude  the  possibility  that  another  detection  method  on  the  subtracted  images  might  improve  the 
detection  capability  of  WASAR. 
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IV.  Findings  and  Conclusions 


4-1  Results  of  This  Research 

The  results  from  this  research  indicate  that  Wide-Angle  Synthetic  Aperture  Radar  definitely 
improves  the  odds  of  detection  over  the  use  of  single  aspect  Synthetic  Aperture  Radar.  In  all  of 
the  detection  methods  used  in  this  research,  pre-set  threshold,  scaled  threshold,  and  subtraction, 
the  use  of  WASAR  always  performed  better  by  detecting  more  of  the  targets  over  a  larger  range  of 
threshold  values. 

One  of  the  discouraging  results  of  this  research,  is  that  the  polarimetric  processing  techniques 
used  in  this  thesis  did  not  seem  to  improve  our  detection  performance  appreciably.  In  most  of 
the  cases  researched,  the  HH  alone  performed  just  as  well  as  any  of  the  polarimetrically  combined 
cases.  This  would  seem  intuitively  incorrect,  since  more  information  is  provided  by  the  three 
polarimetrically  diverse  images,  it  should  improve  the  odds  of  detection.  The  results  did  not 
support  this. 

As  expected,  L-band  (1200  MHz)  radar  performed  well  for  foliage  penetration.  For  most  of 
the  cases  tested  in  this  research,  the  majority  of  the  targets  could  be  detected  even  if  they  were 
hidden  under  foliage.  For  large  amounts  of  foliage,  as  in  scene  three,  the  false  alarm  rate  climbs 
quickly  and  eventually  defeats  our  detection  algorithm  at  the  cross-over  point. 

Lastly,  this  research  indicated  that  subtraction  of  the  aspect  angles  does  show  potential.  The 
subtraction  of  the  aspects  indicates  that  the  targets  do  vary  from  aspect  to  aspect;  while,  the  clutter 
does  not  vary  as  much.  The  results  from  the  pre-set  and  scaled  threshold  detection  methods  also 
support  the  conclusion  that  the  targets  will  vary  with  aspect  angle. 

4.2  Future  Work 

The  results  of  this  thesis  indicate  that  continuing  work  in  this  new  facet  of  Synthetic  Aperture 
Radar  is  very  much  needed.  Numerous  sub-areas  of  research  could  be  explored  more  fully: 
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1.  Use  of  more  sophisticated  detection  techniques  on  the  subtraction  method  discussed  in  the 
previous  chapter,  such  as  the  use  of  a  CFAR  (Constant  False  Alarm  Rate)  detector. 

2.  Use  of  multiple  depression  angle  images.  This  thesis  used  data  sets  that  only  contained  images 
from  a  depression  angle  of  15°.  It  may  be  possible  to  use  the  diversity  between  depression 
angles  in  the  same  manner  as  the  aspect  diversity  was  used  in  this  research. 

3.  Use  of  multiple  frequency  radar  systems.  The  image  sets  generated  for  this  research  used 
only  a  single  frequency,  1200  MHz.  It  may  be  possible  that  probing  a  scene  with  multiple 
frequencies,  with  comparison  of  the  various  returns,  might  yield  some  new  and  better  detection 
techniques. 

4.  Use  of  various  other  target  classes  and  mixture  of  these  classes.  This  thesis  contained  only  one 
type  of  target,  M35  trucks.  With  future  updates  of  the  Loral  software,  more  target  classes  will 
be  available  and  these  new  target  classes  should  be  explored  fully  and  any  existing  detection 
algorithms  should  be  updated  to  include  them. 

5.  Variation  of  seasons  on  the  same  scene.  It  is  possible  that  variations  in  season  will  have  a 
dramatic  effect  on  the  performance  of  any  detection  algorithm.  Since  only  fall  season  was 
explored  in  this  research,  other  seasons  would  need  to  be  explored  as  well  to  verify  any 
detection  algorithm. 

6.  Use  of  netted  and  un-netted  targets  in  the  same  scene.  In  a  real  world  scene,  the  possibility 
definitely  exists  of  having  netted  and  un-netted  targets  within  the  same  area.  Any  detection 
algorithm  developed  would  ideally  be  able  to  pick  out  both  types  of  targets. 

7.  Exploration  of  other  polarimetric  techniques  and  specific  uses  of  HV  and  VV  returns. 

Research  in  this  area  could  also  be  expanded  into  use  of  multiple  sensor  arrays  in  which 
not  only  radar  data  is  used,  but  multi-spectrum,  sonar,  etc.  are  combined  in  an  optimal  target 
detection  algorithm.  Research  in  this  area  may  also  find  application  some  day  in  non-military 
applications  such  as  medical  imaging.  Here  again,  the  list  of  possible  non-military  applications 
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is  endless;  geological  testing,  structural  integrity  testing  for  civil  engineering,  oceanic  biological 
tracking,  etc.  As  can  be  clearly  seen  from  these  incomplete  lists,  this  research  is  definitely  not 
complete  by  any  means,  and  won’t  be  for  many  many  years  to  come. 
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Appendix  A.  Review  of  Loral  CTD  Software 

The  software  that  was  used  in  this  research  to  generate  the  data  was  provide  by  Loral  Corp. 


The  software  was  provided  pre-release  as  a  courtesy.  The  ultimate  purchaser  is  Wright  Laboratory, 
Target  Recognition  Branch.  The  version  that  was  made  available  at  this  time  was  a  stripped  down 
version  of  the  software  that  will  be  made  available  to  Wright  Laboratory  in  the  near  future.  The 
following  paragraphs  describe  some  of  the  capabilities  of  the  software  and  how  it  was  used  for  the 
research  done  in  this  thesis. 

The  SAR  image  data  generation  routine  is  a  user  interactive  VAX- Fortran  program  that 
also  uses  resident  data  files  in  creating  image  data  files.  No  imaging  capabilities  exist  with  this 
particular  software  package.  The  data  is  written  to  complex  VAX-binary  files.  For  our  purposes, 
it  was  necessary  to  convert  this  data  from  binary  into  an  ASCII  file  for  transfer  to  a  Sparc  station 
for  imaging  and  processing  using  the  high-level  language  MATLAB® . 

Once  invoked,  the  data  generation  software  will  query  the  user  for  several  pieces  of  informa¬ 
tion.  The  following  is  a  list  of  the  prompts: 

1.  PV-Wave  available:  This  item  was  not  available  on  the  version  that  was  delivered  to  us. 

2.  Debug  Mode:  The  debug  mode  was  supposed  to  be  available  in  our  version;  however,  though 
it  was  used  a  few  times,  it  did  not  yield  any  information. 

3.  Local  or  regional  statistical  variance:  Only  local  variance  was  available  on  our  version. 

4.  Seed  Value:  This  value  specified  the  seed  value  used  by  the  random  number  generator.  If  this 
value  were  'hanged,  a  different  scene  would  result. 

5-  File  Name:  This  file  name  was  used  as  the  prefix  for  the  data  files  that  would  result.  Each 
polarimetry,  aspect  angle,  depression  angle,  and  frequency  would  have  a  different  suffix  since 
it  was  a  different  image.  The  resultant  output  image  file  contained  a  .img  extension. 
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6.  Polarization  Isolation:  This  is  the  amount  of  separation  that  was  needed  for  the  various 
polarimetries.  A  typical  value  was  -20  dB. 

7.  Percentage  of  foliage  extent:  This  value  would  indicate  how  much  of  the  scene  was  covered  by 
foliage.  A  maximum  of  50%  was  allowed.  If  0%  were  chosen,  the  image  would  be  un-forested. 

8.  Frequency  Mode:  The  number  of  frequencies  that  would  be  used  in  the  simulated  radar 
system.  Choices  were  single,  dual,  triple,  and  quad. 

9.  Frequency  Number:  The  frequency  of  the  probing  radar  system.  Values  could  range  between 
100  to  1500  MHz.  User  would  be  prompted  for  each  frequency  depending  on  number  of 
different  frequencies  chosen  in  last  item. 

10.  Constant  resolution  or  constant  integration  aperture:  This  item  is  fairly  self-explanatory.  For 
the  research  in  this  paper,  constant  resolution  was  chosen  in  order  to  facilitate  registering  the 
images  with  each  other. 

1 1 .  Initial  depression  angle:  This  would  give  the  initial  depression  angle  of  the  airborne  or  elevated 
SAR  system.  Values  of  15°  to  60°  were  allowed. 

12.  Number  of  depression  angles:  The  number  of  different  depression  angles  as  explained  in 
the  previous  item.  The  depression  angles  would  begin  with  the  initial  angle  and  then  be 
incremented  15°  for  each  number  that  is  input  at  this  stage. 

13.  Desired  polarization:  Self-explanatory.  Choices  were  HH,  VV,  HV,  or  fully  polarimetric. 

14.  View  Aspects:  This  would  indicate  the  different  aspect  angles  for  which  a  scene  would  be 
viewed  from.  Choices  were:  single  pass,  +/-  45°,  +/-  30°,  and  +/-  15°.  When  other  than 
single  pass  was  chosen,  the  0°  angle  was  also  generated.  It  was  possible  by  entering  all  other 
parameters  the  same  except  this  one  to  generate  the  same  scene  but  with  more  aspect  angles. 
So  it  was  possible  to  have  up  to  seven  different  viewing  angles. 
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15.  Reference  image  for  change  detection:  This  was  supposed  to  generate  the  same  scene  with  a 
temporal  difference;  however,  on  the  version  that  we  received  the  scene  would  not  change. 

16.  Terrain  type:  Grass  or  short  vegetation  (shrubs).  Only  grass  was  available  for  the  un-forested 
regions  of  the  image. 

17.  TVee  type:  Coniferous,  Deciduous,  or  a  mixture  of  both  could  be  specified. 

18.  Average  leaf  extent:  This  specified  how  much  of  an  individual  tree  would  be  covered  with 
leaves  on  the  average.  A  typical  value  was  7  meters.  Note:  if  0%  foliage  extent  were  chosen 
in  item  7,  this  number  was  ignored. 

19.  Season:  Fall  or  summer  could  be  chosen. 

20.  Large  scene  simulation:  Large  scene  simulation  was  not  available  on  this  version. 

21.  Target  classes:  Allows  the  entering  of  no,  one  type,  or  multiple  types  of  targets.  The  version 
provided  to  AFIT  only  had  available  none  or  one  type.  If  none  was  selected,  then  the  prompts 
would  jump  to  item  (26). 

22.  Target  class:  Targets  were  listed  as  small,  medium,  large,  and  corner  reflectors.  Small  (M35 
trucks)  and  corner  reflectors  were  only  available  with  our  version.  If  the  small  target  class 
is  selected,  the  software  attempts  to  place  up  to  15  targets  in  the  images.  The  number  that 
actually  show  up  in  the  image  is  not  directly  controllable  by  the  operator. 

23.  Netted  targets:  Self-explanatory. 

24.  All  target  conditions  or  deployed  only:  If  deployed  were  selected,  the  targets  would  only 
appear  in  the  forest.  With  all  conditions,  they  may  appear  anywhere  in  the  image. 

25.  Random  or  specified  target  aspects:  Targets  could  only  be  placed  randomly  with  the  version 
available  to  us.  If  targets  were  selected,  then  fifteen  M35  trucks  would  be  randomly  placed 
with  random  aspect  angles. 

26.  Manual  or  random  seeding  of  trees:  Only  random  seeding  of  trees  was  available. 
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27.  Trunk  density  in  trunks  per  square  meter:  Trunk  density  up  to  about  0.01  trunk  per  square 
meter  was  allowed  for  our  version. 

28.  Variable  tree  density:  Self-explanatory. 

29.  Calibration  error:  This  item  not  available. 

30.  ISLR:  Not  available. 

31.  Interference:  Not  available. 

32.  CTD  icon:  Not  available. 

33.  Flight  path:  Not  available. 

34.  Ground  truth:  Would  result  in  a  ground  truth  .pic  image  file.  Though  this  file  was  generated, 
it  could  not  be  displayed  in  the  same  manner  as  the  .img  files;  therefore,  it  could  not  be 
verified  as  to  its  accuracy. 

35.  Output  image  file:  This  would  allow  the  writing  of  the  file  to  the  file  name  specified  previously. 

36.  Integer/Complex/Compress:  Specifies  the  format  of  the  output  file. 

The  software  also  had  safeguards  for  the  most  part  that  would  not  allow  the  input  of  data 
that  was  outside  the  specified  ranges.  However,  when  certain  routines  were  not  available  due  to 
our  having  a  limited  version,  the  routine  would  die. 

The  software  created  several  files  depending  on  options.  The  .lis  files  that  were  created 
contained  a  listing  of  the  target  coordinates  for  that  particular  file.  The  .img  files  were  the  complex 
VAX-binary  image  data  files.  A  .pic  file  may  have  been  created  if  a  ground  truth  image  were 
requested. 

After  the  image  files  were  created,  it  was  necessary  to  put  them  into  some  type  of  format  that 
could  be  transferred  to  the  Sparc  station  for  manipulation  in  M ATLAB® .  This  was  accomplished 
by  an  encoding  routine.  This  routine  would  take  the  binary  files  and  write  them  to  a  .dat  file  that 
was  in  ASCII  format.  This  file  then  could  be  in  theory  read  by  any  machine. 
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The  following  parameters  were  used  for  all  the  data  sets  that  were  generated 

1.  PV-Wave  available:  No 

2.  Debug  Mode:  No 

3.  Local  or  regional  statistical  variance:  Local 

4.  Seed  Value:  57 

5.  Polarization  Isolation:  -20.0 

6.  Frequency  Mode:  single  (1) 

7.  Frequency  Number:  1200 

8.  Constant  resolution  or  constant  integration  aperture:  resolution 

9.  Initial  depression  angle:  15 

10.  Number  of  depression  angles:  1 

11.  Desired  polarization:  fully  polarimetric  (3) 

12.  View  Aspects:  all  were  generated  (2,3,4) 

13.  Reference  image  for  change  detection:  No 

14.  Terrain  type:  grass 

15.  Tree  type:  Deciduous  (2) 

16.  Average  leaf  extent:  7.0 

17.  Season:  Fall 

18.  Large  scene  simulation:  No 

19.  Target  classes:  single  (1) 

20.  Target  class:  small  (3) 

21.  Netted  targets:  No 

22.  All  target  conditions  or  deployed  onl.  all 

23.  Random  or  specified  target  aspects:  random 

24.  Manual  or  random  seeding  of  trees:  random 

25.  Trunk  density  in  trunks  per  square  meter:  0.008 

26.  Variable  tree  density:  yes 

27.  Calibration  error:  no 

28.  ISLR:  no 

29.  Interference:  no 

30.  CTD  icon:  no 

31.  Flight  path:  no 

32.  Ground  truth:  no 

33.  Output  image  file:  no 

34.  Integer/Complex/Compress:  complex 
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The  following  were  the  variant  parameters  for  the  first  data  set: 

1.  File  Name:  tl/ 

2.  Percentage  of  foliage  extent:  0.0 

The  following  were  the  variant  parameters  for  the  second  data  set: 

1.  File  Name:  t2 / 

2.  Percentage  of  foliage  extent:  0.25 

The  following  were  the  variant  parameters  for  the  third  data  set: 

1.  File  Name:  t3/ 

2.  Percentage  of  foliage  extent:  .50 
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Appendix  B.  Computer  Routines 


BA  VAX  File 

The  following  is  a  translation  program  that  took  the  complex  binary  output  of  the  Loral  ctd 
software  on  a  VAX  and  wrote  it  to  an  ASCII  file  for  transfer  to  the  MATLAB®  routines  which 
resided  on  a  Sun  workstation.  This  routine  was  written  originally  by  Steve  Kilberg  of  Loral  and 
modified  by  Patti  Ryan  for  our  specific  needs  at  AFIT. 


C2346678901234687890  Program  XQ_BYT . FOR 

complex  iq(612,612) 

real  max, min 

intsgsr  nl,ns,  lunit 

character *80  inf ils, outfits 

typs  * ,  'Enter  input  filename' 
accept  200,  infile 
200  format (a80) 

type  *,  'Enter  output  filename' 
accept  200,  outfile 

type  ♦,  'Reading  input  file' 
call  read_iq(infile,iq,612) 
type  *,  'Input  file  read’ 

na  =  612 
lout  =  47 

C  open  a  new  file  for  sequential  formatted  access 
C  use  format  statement  300 

open(unit=lout,name=outfile,form='formatted’ ,  status  =  'new’) 

300  format(gl0.3,2x,gl0.3) 

400  format  (/) 


C  another  possibility 
C  300  format(gl0.3) 


do  10  i  =  1,612 

uritedout,  300)  (iq(j,i),  j  =  1,612) 
write(lout,400) 

10  continue 


C  try  opening  a  new  file  for  unformated  access 
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C  open(unit*lunit ,  nauesoutfile,  for*=’uniormatted' ,  status  =  'nee') 
C 

C  print  oat  loop  lor  unformatted 

C  do  10  i  *  1,512 
C  srite(6,*)(iq(j.i),  j=l,512) 

C  10  continue 

close  (limit) 

stop 

end 

c************************************************************ 

subroutine  read_iq(inlile, array, iugsize) 

coupler  array (612, 5 12) 

integer  luin.nline , nsanp , iugsize 

character*80  infile 

nline  =  ingsize 
nsanp  =  iugsize 


lain  =  66 

open  (unit  slain ,  naue= inf  ile ,  fom=  *  unformatted 1 , 
ft  access* 'direct' .recordtype* 'fixed* , status  *  'old') 

do  i  *  1, nline 

read(luin  ’i)  (array(j,i),  j  *  l.nsaup) 
end  do 

close (unit  *  lain) 

return 

end 

<>********«************************************************** 
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B.t  Sun  Files 


The  scripts  contained  in  this  appendix  are  the  MATLAB®  routines  that  were  used  to  load, 
process,  and  display  the  image  data  that  was  generated  by  the  Loral  software.  The  first  routine 
is  what  could  be  called  the  main  program.  This  is  the  routine  that  is  invoked  every  time,  sar.m 
contains  the  menu  that  calls  all  the  other  major  routines.  Each  file  contains  a  header  which 
explains  the  purpose  of  the  file.  Comments  on  particular  phases  of  the  code  were  not  added  where 
text  output  messages  explained  the  procedure.  Maximum  run  times  are  hard-coded  and  would 
necessarily  be  cpu  load  dependent.  Quirks  do  not  necessarily  indicate  problems  with  the  code,  they 
might  simply  be  items  of  particular  interest  for  that  routine  or  function. 

MATLAB®  does  have  an  imaging  quirk  that  was  discovered  during  this  research.  MATLAB® 
does  not  display  images  in  the  same  coordinate  system  as  it  plots.  When  using  the  “image()” 
function  to  display  the  image  matrices,  the  displayed  figure  will  transpose  the  image.  But,  when 
a  plot  of  target  coordinates  is  overlayed  on  this  image,  the  coordinates  are  plotted  true  to  axis 
coordinates.  So,  a  plot  and  an  image  in  the  same  figure  will  not  use  the  same  coordinates.  In 
order  to  correct  for  this  idiosyncrasy,  it  was  necessary  to  transpose  the  image  before  displaying  so 
that  plotting  could  be  overlayed  on  the  images.  For  the  particular  coding  involved,  observe  the  file 
called  “display .m”. 

The  code  that  has  been  written  for  this  research  is  not  optimal  in  any  sense  of  the  word. 
The  code  was  written  for  my  understanding  primarily,  for  optimization  of  available  memory  (swap 
and  disk)  secondly,  and  speed  lastly.  Since  the  purpose  of  this  thesis  was  conceptual,  coding 
optimization  was  not  given  a  high  priority;  however,  since  MATLAB®  is  a  high-level  C  language, 
direct  conversion  should  be  no  problem.  After  conversion,  optimization  might  best  be  handled  by 
a  software  engineer. 


B-3 


%•*********************************•****************************************** 
X  RDUTIIE  ur.i 

%*************«****************************•********************************** 
X  Purpose:  Main  aenu  rout  in® .  Create®  a  figure  for  display  of  aenu  and  other 
X  data. 

X 

X  Called  by:  KatLab 

X  Calls:  filesl.a,  input.data.a,  filt.a,  detection. a,  show, s tat. a, 

X  displayl.a,  display2.a,  dell. a,  del2.a 

X 

X***************************************************************************** 
clear  all 

X  check  for  existence  of  "files. aatH;  if  it  doesn't  exist,  create  it 
if  exist ('files.aat')  "=  2 
files 1 

save  files  directory  file_naae  prep  detect 
end 
clear 

X  Menu  title 

f igureC 'Position* , [635.5,516,400]  .... 

'lane' , 'SAR  Processing' , 'VuaberTitle' , 'off ’) ; 


X  aenus 

aenul  =  uiaenu ( 'Label' , '  Input  '); 

uiaenu(aenul, 'Label' , 'Input  .dat  files' , 'Callback' , 'input _data') ; 

aenul  *  uiaenu ( ' Label ’ , '  Pre-processing  '); 

uiaenu (aenul , 'Label' , 'lone' , 'Callback' , 'filt(O) ' ) ; 

uiaenu  (aenul,  'Label' ,  'Rotation  Only' ,  'Callback' ,  ’filt(l) ') ; 

uiaenu (aenul, 'Label' , ’Average’ , 'Callback' , ’filt(2) ') ; 

uiaenu(aenul ,  'Label' ,  'Span' , ' Callback* ,  ’iilt(3) ' ) ; 

uiaenu  (aenul,  'Label' ,  'Optiaal  Weighting',  'Callback' ,  *filt(4)  ’); 

uiaenu  (aenul ,  'Label' ,  'Polariaetric  Whitening' ,  'Callback' ,  'filt(S) ') ; 

uiaenu (aenul , 'Label ’ , ’ Create  all ' ,  ... 

'Callback' , 'f ilt(0) ;filt(l) ;filt(2) ;filt(3) ;f ilt(4) ;filt(5) ’) ; 

aenul  -  uiaenu ( 'Label* , ’  Detection  '); 

uiaenu (aenu 1 , 'Label' , 'lone' , 'Callback' , ’detection(O) '); 

uiaenu  (aenul ,  'Label* ,  'Multiple  Detection  (preset  threshold)’, _ 

'Callback' , 'clear  all;detection(2) ’); 

uiaenu(aenul, 'Label' , 'Multiple  Detect  (scaled  threshold)  ’ _ _ 

'Callback',  ’dear  all;d«u«ction(l) ’) ; 
uiaenu(menul , 'Label ’ , 'Subtraction' ,  ... 

'Callback' ,  'dear  all;detection(3)  ’); 

aenul  *  uiaenu ( 'Label' , *  Statistics  '); 

uiaenu  (aenul , 'Label* , 'Statistics  on  inages’ , 'Callback' , 'shou.stat') ; 
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■anal  ■  uimenu ('Label' , '  Output  ’); 

uimenu(menul , 'Label ’ , ’Display  images' , ’Callback' , 'display 1(0) ’) ; 
uimenu(menul , ’Label’ , ’Display  prf  images' , ’Callback' , 'displayl(-l) ' ) ; 
uimenuCmenul , ’Label ' , 'Display  images  (autoscale) ' , ’Callback’ , 'display 1(1) ' ) ; 
uimenu (menul, 'Label’ , ’Create  eps  file’ , 'Callback’ , ’display2') ; 

menul  =  uimenu( ’Label' , ’  Exit  ’); 

uimenu(menul, ’Label' ,  ’Exit  from  Natlab', 'Callback', 'quit'); 

uimenu (menul, ’Label* , ’Exit  to  Natlab* , ’Callback’ , ’dell; clear ;delete(l) ') ; 

uimenu (menul , 'Label ' , ’Reset ’ ,  ... 

'Callback* , *dell;clear;delete(l) ;sar; ’) ; 

menul  =  uimenu ( ’Label* , ’  Cleanup  *); 

uimenu(menul ,  ’Label’ ,  ’Clear  figures* ,  ’Callback' ,  ’dell*) ; 
uimenu (menul, ’Label* , *  Delete  .mat  files' , ’Callback’ , *del2; *) ; 

clear  menul 

%»**************************************************************************** 


%*************************************«***************♦*************«********* 
X  ROUTINE  filesl.m 

X***************************************************************************** 
X  Purpose:  Used  to  initialize  files. mat.  This  file  must  be  changed  uhen  a 
X  new  data  set  is  being  loaded.  All  .mat  files  need  to  be  deleted 

X  from  current  directory. 

X 

X  Called  by:  sar.m 
X 

X  Quirks : 

X  1)  Variables  are  hard-coded.  File  must  be  edited  to  update 
X 

X***************************************************************************** 

directory  =  ’/tap_ant/hoae/hawkeye8/kknurr/images/ctdimages’ ; 
file.name  *  C’t2_p45’  ;  *t2_p30*  ;  *t2_pl5’  ;'t2_000’  ;  ... 

*t2_ml5*  ;  *t2_n30*  ;  ’t2_m45'J; 

prep  *  ’non’ ; 
detect  =  0; 

%***************************************************************************** 
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X***************************************************************************** 
X  ROUTIIE  inpat .data. ■ 

X***********************************************************e***************** 
X  Pnrpoaa:  Loads  data  Iron  . dat  (ASCII)  files  and  creates  an  annotated 
X  inage  matrix.  Also,  loads  target  coordinates.  Fully  polar 

X  inages  are  then  stored  in  one  of  seven  unrot at edX. mat  files 

X 

X  Called  by:  sar.m 
X  Calls:  targlis.n 
X 

X  Quirks : 

X  1)  Directory  changes  are  hard-coded 
X 

%****************************************************,t*******4tt,t4„^,4t4t** 

tine  =  f  ix(dock)  ; 

X  setup  message  window 
figure(l) 
cla 

hold  on 

axis('xy’ , 'off’) 
nessl  =  text( . 1, .6, ’  ’); 
tness  =  text (.1,1, '  ’); 
aness  =  text( . 1, .95, ’  '); 

set  (tness,  'String'  .[’Starting  Tine:  ’  nun2str  (tine  (1,4))  ... 

nun2str(tine(l ,5))  ’:'  nun2str(tine(l,6))  ]); 
set (aness, 'String', 'Max  Run  Tine:  90  minutes');  pause(l) 

if  exist (’files. mat')  '=  2 

set(nes8l, 'String' , 'Files. mat  information  must  be  loaded  first') 
pause (5) 
figured) 
cla 

hold  off 
else 

load  files 
prep  =  'non'; 

save  files  file.name  directory  prep  detect 

X  load  iaages  three  at  a  time  (one  aspect  at  a  time) 
for  i  =  1  :  sized ile.name,  1) 

set (nessl, 'String', [’Loading  ’  file_name(i, :)  '  Target  Coordinates’]) 
pause(2) 

targets  =  targlis(  [file_nane(i, :)  ’.dat.lis’]  ); 

X***********change  directories  (hard  code)********** 
cd  / tmp.mnt /home/havkey e8/kknurr /inages/ctdimages 
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X*************************************************** 

set (ness 1,  'String*.  ['Loading  *  f ile_nane(i, :)  *  EH  polar  iaetry’]  ) 
pause(l) 

fid  *  fopen([f ile_naae(i, :)  ’_hh.dat’] , ’r’ ) ; 
hh  *  fscanf (fid, 'Xf ’); 
status  =  fcloss(fid) ; 

s*t(a*ssl, ’String* , ['Loading  ’  file_nane(i, : )  ’  HV  polarinetry']  ) 
panss(l) 

fid  =  f op«n( [f ile_naae(i , : )  ’_hv.dat’] , ’r’) ; 
hv  =  fscanf(fid.’Xf’); 
status  =  fcloss(fid); 

sat (ass si,’ String’, [’Loading  ’  file_naae(i, :)  ’  W  polariastry’]  ) 
pauaa(l) 

fid  =  fopan([fila_naaa(i, :)  ’_w.dat’]  ,’r’); 
nr  =  fscanf (fid, ’Xf '); 

status  =  f clo8a(f id) ; 
claar  fid  status 

X***********change  diractorias  (hard  coda)********** 
cd  /tap_ant/hoas/haTkaye8/kknurr/iaagas 
X*** ***********♦*♦*♦♦♦♦♦♦******♦♦♦♦♦♦♦♦♦*♦♦♦♦ *♦♦♦♦*# 

sat(aassl, ’String’ .['Creating  ’  file_nane(i, : )  ’  HH  iaage  file’]  ) 
pausa(l) 

hh  =  hh( 1:2: langth(hh) -1)  +  j*hh(2:2:length(hh)) ; 

hh  *  reshape(hh,sqrt(length(hh)).sqrt(length(hh))); 
sat (aesal,’ String’, f Creating  •  iile_nune(i, :)  ’  HV  inage  file’]  ) 
paused  ) 

hv  =  hv(l: 2: length (hv)-l)  +  j*hv(2:2:length(hv)) ; 
hv  *  reshape  (hv ,  sqrt  (length  (hv)  ) ,  sqrt  (length  (hv)  )  ) ; 
set(aessl, ’String’, [’Creating  ’  file_naae(i, : )  ’  W  iaage  file’]  ) 
paused) 

w  =  vv(l:2:length(w)-l)  +  j*vv(2:2:langth(vv)) ; 

W  =  reshape(vv, sqrt (length(vv) ) , sqrt (length (vv)))  ; 

set(aessl, ’String’, [’Saving  as  unrotated’  nua2str(i)  ’.mat’]  ) 
paused) 

X  Select  file  name 
if  i  ==  1 

save  unxotatedl  hh  hv  w  targets 
elseif  i  ==  2 

save  unrotated2  hh  hv  vv  targets 
elseif  i  ==  3 

save  unrotated3  hh  hv  w  targets 
elseif  i  ==  4 

save  unrotated4  hh  hv  vv  targets 
elseif  i  ==  S 

save  unrotated6  hh  hv  vv  targets 
elseif  i  ==  6 

save  unrotatedS  hh  hv  vv  targets 
elseif  i  ==  7 

save  unrotated7  hh  hv  w  targets 
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and 

dear  hh  hr  vv  targets 
and 

set(aessl , 'String' , 'Input  of  files  complete ’) 
paase(6) 
figured) 
cla 

hold  off 
end 
clear 

%******************************************************************* 
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X****************************************************************************** 
X  FUICTIO*  targlis.a 

X****************************************************************************** 

X  Purpose:  Loads  data  from  .lis  (ASCII)  files  and  creates  an  unrotated 
X  target  coordinate  matrix.  Strips  off  header  and  retains  only  the 

X  coordinate  information  from  file  as  it  is  formated  by  ctd 

X  software. 

X 

X  Called  by:  input _data.m 
X 

X  Quirks : 

X  1)  Directory  changes  are  hard-coded 

X*************«>*************************************************************** 

function  [targets]  =  targlis(file.name) 

X***********change  directories  (hard  code)********** 
cd  /tmp_mnt /home/hawkey e8/kknurr/images/ctdimages 

x*************************************************** 

fid  =  fopen(f ile_name , 'r ‘ ) ; 

X  strip  off  header 
read_in  -  *n*; 
while  read_in  '=  • ; ' 

read_in  =  fscanf (f id, ’Xs* ,1) ; 
header  =  [header  read. in] ; 
end 

read_in  =  fscanf (fid, 'Xs' ,6) ; 

X  read  in  target  coordinates 

read_in  =  fscanf (fid, 'Xs' , 1) ; 
while  read.in  "=  □ 

targets  =  [targets  ;  fscanf (fid, 'Xf' ,1)  fscanf (fid, ’Xf’ , 1)]  ; 
read.in  =  fscanf (fid, ’Xs’ ,10) ; 
read_in  =  fscanf (fid, ’Xs’ ,1); 
end 

status  =  fdose(fid); 

X***********change  directories  (hard  code)********** 
cd  /tmp_mnt/home/hawkeye8/kknurr/iaages 
X*************************************************** 

clear  fid  status  header  read.in 

X****************************************************************************** 
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%***************************************************************************** 
X  BOOTHE  rot.  a 

x***************************************************************************** 

X  Purpose:  Rotates  .nat  files  that  are  created  by  input  data. a. 

X 

X  Called  by:  sar.si 
X  Calls:  rotation. a 
X 

X  quirks : 

X  1)  Rotation  aaount  is  hard-coded  for  7  aspect  angles  at  IS  degree 
X  increuents.  iutoaatic  coding  would  not  work  (see  anotated  code). 

X 

X************************************************ ******** ********************* 

tine  =  fix(clock); 

X  setup  nessage  windows 
figured) 
cla 

hold  on 

axis(’xy’ ,  'off') 
aessl  =  text(.l,.6,’  ’); 
taess  =  text(.l,l,’  ’); 
aaess  =  text(. 1, .95, ’  ’); 
trmess  »  text(.l,.9,’  ’); 
armeas  *  text(.l, .86,’  ’); 

set  (taess,  ’String' ,  ['Starting  Tine:  ’  nua2str(tiae(l,4))  ... 

nua2str(tiae(l ,5))  ':’  nua2str(tiae(l ,6))  ]); 
set (aaess, ’String’ , ’Max  Run  Tine:  260  ainutes'); 
set (araess, ’String' , ’Max  Rotation  Tine:  13  ainutes');  pause(l) 

load  files 

if  exist (’unrotatedl. aat’)  '=  2  X  else  go  to  end 

8et(aessl,  'String' ,  'Unrotated  files  are  not  in  this  directory’) 
pause (5) 
figure(l) 
cla 

hold  off 
else 

X  load  iaages  three  at  a  time 
for  i  =  1  :  sized ile_name,  1) 

set (aessl, 'String' ,  ['Loading  unrotated’  num2str(i)  '.mat']  ) 
paused ) 
if  i  ==  1 
load  unrotatedl 
elseif  i  ==  2 
load  unrotated2 
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elseif  i  ==  3 
load  unrotated3 
elseif  i  ==  4 
load  unrotated4 
elseif  i  «  5 
load  unrot at  ad  6 
elseif  i  ==  6 
load  unrotated6 
•laai t  i  ==  7 
load  unrotated7 
and 

X******************************************************************** 

X  Tha  following  rotation  angla  determination  is  broken,  angle  must  be 
X  hardcoded  at  this  point.  This  code  attempted  to  use  the  filename 
X  information  to  determine  amount  of  rotation. 
%******************************************************************** 

X  Determine  angle  of  rotation 
X  if  i  >  1 

X  for  j=l:size(file_name,2) 

X  if  file_name(i, j)==’_’  ... 

X  *  (  file_name(i, j+l)==’p’  I  f ile_name(i, j+l)==’0’  ) 

X  number  =  abs(f ile_name(i, j+2: j+3)); 

X  number  ■  round (abs(  abs(file_naae(l, j+2: j+3))  -  number  )); 

X  elseif  f ile_name ( i , j ) == ’ _ ’  *  file_name(i, j+l)=='m’ 

X  number  =  -abs(file_name(i, j+2: j+3)); 

X  number  =  round (abs(  abs(file.name(l, j+2: j+3))  -  number  )); 

X  end 

X  end 

X  else 
X  number  =  0; 

X  end 
Xif  i  ==  1 
X  number  =  0; 

Xelseif  i  ==  2 
X  number  =  3; 

Xelseif  i  ==  3 
X  number  =  6; 

Xend 

number  =  i-i;  X  hard-coding  for  7  aspect  angles 

%**4  j************************************************** 

X  rotate  over  "number"  of  15  degree  increments 
time  =  fix(clock); 

set  (trmess,  ’String' ,  [’  Rotation  start  Time:  ’  num2str(time(l,4))  ... 

num2str(time(l,5))  num2str(time(l,6))  ],  ... 
’Visible’ , ’on’) ; 

set(messl, ’String’ , ['Rotating  ’  file_name(i, :)  ’  target  coordinates']) 


B-ll 


paus«(l) 

teap  =  0  *  ones (512, 512); 
for  j=l:size(targets,l) 
if  targats(j.l)  *=  0 
tmpC  targsts(j,l),targets(j,2)  )  =  255; 
rad 
and 

trap  =  rot at ion (trap, number ) ; 

[k.l]  =  f ind(trap==255) ; 
targets  =  [  k  1  ]; 
clear  j  k  1  trap 

tiae  =  f  ix(dock) ; 

set  (traess,  ’Stris^*,  ['Rotation  start  Tiae:  ’  nua2str(tiae(l  ,4)) 
uua2str(tiae(l,5))  *:*  num2str  (timed ,6))  ]); 
aet(messl, ’String’ , [’Rotating  ’  file_naae(i, : )  *  HH  image  file’]) 
pauae(l) 

hh  =  rotation(hh, number) ; 
tiae  =  fiz(dock); 

set  (traess,  ’String’ ,  [’Rotation  start  Tiae:  ’  nua2str(tiae(l,4))  ’:’ 

nua2atr( tiae(l ,5))  ’:’  nua2str( timed ,6))  ]); 
set(aessl, ’String’ , [’Rotating  ’  f ile_naae(i, : )  ’  HV  image  file’]) 
pause(l) 

hv  =  rotation(hv, number) ; 
time  =  f iz(clock) ; 

set (traess, ’String' , ['Rotation  start  Time:  ’  nun2str(time(l,4)) 
nua2str(time(l ,5))  ’:’  nua2str(time(l ,6))  ]); 
set(aes8l, ’String’ , [’Rotating  ’  file_naae(i, : )  ’  W  image  file’]  ) 
pause(l) 

▼v  =  rot at ion (vv, number) ; 
set (traess , ' Visible ’ , ’ of f ’ ) ; 


set (aessl, ’String’ , [’Saving  as  rotated’  nua2str(i) 
pause(l) 


if  i  ==  1 


save 

elseif 

save 

elseif 

save 

elseif 

save 

elseif 

save 

elseif 

save 

elseif 

save 


rotatedl  hh  hv 
i  ==  2 

rotated2  hh  hv 

i  ==  3 

rotated3  hh  hv 
i  ==  4 

rotated4  hh  hv 
i  ==  5 

rotatedS  hh  hv 
i  ==  6 

rotated6  hh  hv 
i  ==  7 

rotated7  hh  hv 


vv  targets 
vv  targets 
vv  targets 
vv  targets 
vv  targets 
vv  targets 
vv  targets 


’.nat’]  ) 
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•ad 

clear  hh  hr  n  target* 

•ad 

■et(aee*l, 'String' , 'Rotation  of  file*  complete' ) 
pau*e(5) 
ligure(l) 
da 

hold  oil 
prop  =  'rot' ; 

save  files  directory  lile_naae  prep  detect 
end 

%*********************************************************«***********«******* 
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%«**************************************************************************** 
X  FU1CTI0*  rotation. ■ 

%*************************************e*********************************«***** 
X  Pnrpoao:  Rotates  an  image  input  by  "number"  of  fifteen  degree  increments. 

X  Also  cuts  image  down  to  overlapping  square  region. 

X 

X  Called  by:  rot. a 
X 

X  Quirks : 

X  l)  Uses  short  cut  for  0  and  +/-  90  degree  rotations. 

X 

X***************************************************************************** 

function  [temp2]  =  rotation(templ .number) 

X  number  indicates  the  number  of  IS  degree  increment  rotations, 
phi  *  number  *  (-pi/12);  X  angle  of  rotation 

X  (a.b)  central  point  of  rotation 
a  =  round  (  size  (tempi,  0/2); 
b  =  round(size(teapl,2)/2) ; 

if  phi  ==  0 

temp 2  =  tempi;  X  0  degree  short  cut 
elseif  phi  ==  pi/2 

temp2  =  rot90(templ  ,-l) ;  X  90  degree  short  cut 
elseif  phi  ==  -pi/2 

temp 2  =  rot90 (tempi, 1) ;  X  90  degree  short  cut 
else 

temp2  =  min(min(templ))  *  ones(size(templ,l),size(templ,2)); 
for  i=l:size(temp2,l) 
for  j=l:size(temp2,2) 

r  =  sqrt((i-a)*2  +  (j-b)*2); 
if  i==a  k  j>=b 
theta  =  0; 
elseif  i*=a  *  j<b 
theta  =  pi; 
else 

theta  =  atan2((j-b),(i-a)); 
end 

x  =  round (r  *  cos (theta  +  phi))  +  a; 
y  =  round (r  *  sin(theta  +  phi))  +  b; 
if  x>=l  t  x<=size(templ,l)  k  y>=l  k  y<=size(templ,2) 
temp2(i, j)  =  templ(x.y); 
end 
end 
end 
end 

X  cutting  down  image  to  square  image  of  overlapping  region 
1  =  a  -  round(sqrt(  (a“2)  /  2  )); 
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a  •  •ix*(t««pl,l)  -  1; 
t«ap2  *  t«ap2(l:a,l:u) ; 

%**#**********************************************************v*************** 


%***************************************************************************** 
X  FUICTIOI  flit.* 

x***************************************************************************** 

X  Purpose:  Selects  the  type  of  pre-processing  that  ia  aaad.  If  the 
X  pre-proceaaed  filea  do  not  exist,  it  craataa  the*. 

X 

X  Called  by:  ear.* 

X  Calls:  ave.a,  span.*,  ous.a,  psf.* 

X 

X***************************************************************************** 

function  □  =  f ilt(neuprep) 

tiaa  =  fix (clock) ; 

X  setup  aessage  uindou 
figured) 

da 

hold  on 

axis('xy' ,  'off') 
see a  =  text (.1,. 6,’  *); 
taess  =  text(.l,l,'  ’); 
aaesa  =  text( . 1, .95, ’  '); 

setCtaesa,  ’String’ ,  [’Starting  Ti*e:  *  nua2str(ti*e(l,4))  ... 

nua2str (time (1,6))  nu*2str (time (1,6))  ]); 

set (aaesa, 'String* ,  'Max  Run  Tiae:  5  ainutea');  pause(l) 

load  files 

if  neuprep  ==  0 
prep  =  ’non* ; 

if  exist Cunrotatedl.aat')  *=  2 
input.data 
end 

elseif  neuprep  ==  1 
prep  =  'rot* ; 

if  exist (’rotatedl.aat*)  *=  2 
rot 
end 

elseif  abs (neuprep)  ==  2 
prep  =  ’axe* ; 

if  exist ('averaged. mat’)  *=  2 
neuprep  =  -1; 
end 

elseif  abs (neuprep)  ==  3 
prep  =  'spn'; 

if  exist ( 'spaned.Bat')  ~=  2 
neuprep  =  -1; 
end 

elseif  abs (neuprep)  “  4 
prep  =  'ous’ ; 
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if  exist (’ opt iaal.aat')  "=  2 
nexprep  *  -1; 
end 

elseif  abs (nexprep)  5 

prap  =  ’pxf  ’ ; 

if  existOehiteaed.aat*)  '=  2 
nexprep  =  -1; 
and 
and 

sax*  files  directory  file.naae  prep  detect 

if  nexprep  <0  X  process  the  files,  otherwise  files  eill  be  loaded 

if  prep==’ave’  I  prep== ’ spn ’  |  prep=='osa’  I  prep=='pwf' 
for  i  =  l:size(file_naae,l) 

set  (mess. 'String* Loading  rotated  *  file_naae(i. :)]) 
panse(l) 
if  i  ==  1 
load  rotatedl 
elseif  i  ==  2 
load  rotated2 
elseif  i  ==  3 
load  rotated3 
elseif  i  ==  4 
load  rotatedd 
elseif  i  ==  5 
load  rotatedS 
elseif  i  ==  6 
load  rotatedfi 
elseif  i  *=  7 
load  rotated7 
end 

set(aess, 'String’ , [’Pre-processing  •  file_naae(i, :)]) 
paused  ) 
if  prep  ==  ’axe' 

tasrp  =  axe(hh.hx.xx) ; 
elseif  prep  ==  ’spn* 
teap  =  span(hh.hx.xx) ; 
elseif  prep  ==  ’oxs' 
teap  =  oxs(hh,hx,xx); 
elseif  prep  ==  ‘pvt’ 
teap  =  pwf (hh.hx.xx) ; 
end 

if  i  ==  1 

aspectl  =  teap; 
elseif  i  ==  2 
aspect 2  =  teap; 
elseif  i  ==  3 
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aspect3  =  trap; 
elseif  i  ==  4 
aspect4  s  trap; 
elseif  i  =3  S 
aspects  3  trap; 
elseif  i  33  e 
aspects  3  teap; 
elseif  i  33  7 
aspaetT  *  trap; 
end 

clear  trap 

and 

and 

if  prap  ==  ’ave’ 

sat(aass, ’String' , ’Saving  as  averaged') 
paused  ) 

save  averaged  aspectl  aspect2  aspect3  aspect4  . . . 

aspects  aspects  aspect?  targets 
elseif  prep  33  >spn’ 

set (aess, 'String* , ’Saving  as  spaned') 
paused) 

save  spaned  aspectl  aspect2  aspect3  aspects  . . . 

aspects  aspects  aspect7  targets 
elseif  prep  ==  'ows ’ 

setCaess , ’String* , ’Saving  as  optimal’) 
paused  ) 

save  optimal  aspectl  aspect2  aspect3  aspect4  . . . 

aspects  aspects  aspect7  targets 
elseif  prep  ==  'pwf ’ 

set (ness , ’ String ’ , ’ Saving  as  whitened ’ ) 
panse(l) 

save  whitened  aspectl  aspect2  &spect3  aspects  . . . 

aspects  aspects  aspect7  targets 

end 

end 

X  Finish  up 

set (mess, ’String' , ’Processing  Complete’) 
pause (6) 
figure(l) 
da 

hold  off 

%***************************************************************************** 
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%************«*************«*************************************************** 

X  FUICTIO*  tTt.i 

x****************************************************************************** 

X  Purpose:  Averages  the  polar imetric  images  to  create  a  nee  image. 

X 

X  Called  by:  filt.m 
X 

X***************************************************************************** 
function  [average]  =  ave(hh,hv,vv) 
average  =  (hh  +  hv  +  vv)/3; 

y,****************************************************************************** 


%*********************************************«******************************** 

X  FUICTIO*  apan.m 

%******+*********************************************************************** 

X  Purpose:  Createa  a  span  image  from  the  polarimetric  images. 

X 

X  Called  by:  filt.m 
X 

%********e******************************************************************** 
function  [sp]  =  span(hh,hv,vv) 

ap  =  abs(hh)."2  +  2*(abs(hv) .*2)  +  abs(w).“2; 

%****************************************************************************** 
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%************************************«***************************************** 
%  FU1CTI01  oss.m 

%***********************************************♦****************************** 
X  Purpose:  Creates  an  Optimally  Weighted  Sun  of  the  polar  inages. 

X 

X  Called  by:  filt.n 
X 

X***************************************************************************** 

function  [y]  =  ous(hh,hv,vv) 

rho  =  nean(mean(hh. *conj (vv) ) )  ... 

/sqrt(  nean(nean(abs(hh) . *2))  .*  nean(nean(abs(vv) . ~2))  ); 
epsilon  *  mean(mean(abs(hv) . ~2))  /  nean(nean(abs(hh) . *2)) ; 
gamut  =  nean(nean(abs(vv) . ~2) )  /  nean(nean(abs(hh) . "2)) ; 
k2  *  (1  +  abs(rho)*2)/epsilon; 
k3  =  1/ganna; 

y  =  abs(hh).~2  +  k2*(abs(hv) . *2)  +  k3*(abs(vv) . “2) ; 
X*************************************’*************'*************'*************** 


X****************************************************************************** 
X  FUICTIO*  pvf . n 

X****************************************************************************** 
X  Purpose:  Creates  a  Polarinetrically  Whitened  inage  fron  the  polar  inages 
X 

X  Called  by:  filt.n 
X 

X***************************************1************************************** 

function  [y]  =  pwf (hh.hv, vv) 

rho  =  nean(nean(hh. *conj (vv) ) ) . . . 

/sqrt(  nean(nean(abs(hh).“2))  .*  mean(nean(abs(vv) . *2))  ); 
epsilon  =  mean(nean(abs(hv) . ~2))  /  mean(nean(abs(hh) . *2)) ; 
ganna  =  mean(nean(abs(vv) .‘2))  /  mean(nean(abs(hh) . “2)) ; 
sigma  =  mean(nean((abs(hh)) . *2)) ; 

kl  =  1/ (sigma  *  (1  -  (abs(rho))‘2)) ; 

k2  =  1/ (sigma  *  epsilon); 

k3  =  1/ (sigma  *  (1  -  (abs(rho))*2)  *  gamma); 

k4  =  (2  *  abs (rho))/ (sigma  *  (1  -  (abs(rho))“2)  *  sqrt(gamma) ) ; 

y  =  kl*(abs(hh)).‘2  +  k2*(abs(hv)) . *2  +  k3*(abs(vv)) . ~2  ... 

-  k4*(abs(hh).*abs(w) .*cos(angle(hh)-angle(vv)-angle(rho))) ; 


X****************************************************************************** 
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X***************************************************************************** 
X  FUICTIOI  display l.n 

%***************************************************************************** 
X  Purpose:  Displays  the  current  pre-processed  inage  data  set  on  the  screen. 

X  Also  sill  indicate  known  target  coordinates  in  selected  images. 

X 

X  Called  by:  sar.a 
X  Calls:  display. n,  dell.n 
X 

X  Quirks : 

X  1)  If  data  set  uses  more  than  three  aspects,  figures  sill  overlap 
X 

%*************•*************************************************************** 

function  □  =  display 1 (auto) 

X  clearing  figures 
dell 

X  setup  message  window 
figure(l) 

da 

hold  on 

axis('xy’ , ’off’) 

mess  =  text (.1,  .6, ’broke  before  anything  ran’); 
load  files 

X  choose  between  polarinatric  and  processed  images 

if  (prep==’non’  I  prep==’rot’)  k  detect* =3  X  polarinetric 

for  i  =  l:size(file_name,l) 

if  prep  ==  'non* 

set  (mess,  ’String’ ,  [’Loading  unrotated  ’  file_name(i, :)]) 
pause(l) 
if  i  ==  1 

load  unrot at edl 
elseif  i  ==  2 
load  unrotated2 
elseif  i  ==  3 
load  unrotated3 
elseif  i  ==  4 
load  unrotated4 
elseif  i  ==  5 
load  unrot at edS 
elseif  i  ==  6 
load  unrot at ed6 
elseif  i  ==  7 
load  unrot at ed7 
end 
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elseif  prep  =5  'rot’ 

set (mss,  ' String* .  ['Loading  rotated  *  file_naae(i, : )] ) 
paused) 
if  i  ==  1 
load  rotatsdl 
elseif  i  ==  2 
load  rotat«d2 
elseif  i  3 
load  rotatod3 
elseif  i  ==  4 
load  rotated4 
elseif  i  ==  5 
load  rotatedS 
elseif  i  as  e 
load  rotatedS 
elseif  i  as  7 
load  rotated7 
end 
end 

if  size(file_naae, 1)  >  3 

rect  =  [(i-l)*105+6, 486,200, 200]; 

else 

rect  =  [(i-l)*210+5,485, 200,200]; 
end 

set(aess , ’String* ,  ['Displaying  *  file_naae(i, :)  ■  hh']) 
paused ) 

hand  =  figure( 'Position' , rect,  ... 

’lane’ . Cf ile_name(i, :)  *  hh’] , 'InaberTitle' , 'off ’) ; 
display (hh , hand . auto) 

X  Indicate  targets 

for  j=l :length(targets) 
if  targets(j,l)  '=  0 

text (targets (j ,1) , targets (j ,2) , 'X') 
end 
end 

if  size(file_naae,l)  >  3 
rect  =  [(i-l)*106+5, 485-230, 200, 200] ; 

else 

rect  =  [(i-l)*210+5, 485-230, 200, 200] ; 
end 

set(aes8, 'String' , [’Displaying  ’  file_naae(i, :)  ’  hv’]) 
pause(l) 

hand  =  figure ( 'Position' , rect,  ... 

’■a»e' , [file_naae(i, :)  ’  hv'] , ’InaberTitle’ , 'off ’) ; 
display (hv , hand , auto) 

if  size(file_naae,l)  >  3 
rect  =  C(i-l)*105+5, 466-460, 200, 200]; 

else 
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rect  =  [(i-l)*210+S, 465-460, 200, 200] ; 
end 

set (■•ss, ’String* , ['Displaying  ’  file_name(i, : )  '  vv’]) 
pan*«(l) 

hand  =  figure( ‘Position' ,rect, 

’lame ’ , [f ile_name(i , : )  *  vv']  ,  'lumberTitle' ,  ’off’) ; 
display ( vv , hand , auto) 

end 

elseif  prep== ’ ave ’  I  prep==’spn'  I  prep== ’ ows '  |  prep=='pwf ’  I  detect 

if  detect  "=  3 
if  prep  ==  ’ave’ 

set  (ness ,  ’String’ ,  ’Loading  averaged’ ) 
pause(l) 
load  averaged 
elseif  prep  ==  ’spn’ 

set(mess, ’String’ , ’Loading  spaned’) 
pause(l) 
load  spaned 
elseif  prep  ==  ’ows’ 

8et(ness, ’String' , 'Loading  optisial’) 
pause(l) 
load  optinal 
elseif  prep  ==  ’pwf ' 

set (mess, 'String' , 'Loading  whitened*) 
pause(l) 
load  whitened 
end 

elseif  detect  ==  3 

setCmess, ’String’ , 'Loading  Subtracted') 
pause(l) 
load  det3 
end 

for  i  =  1 :size(file_naae,l) 
setCmess,  ... 

’String’ , [’Displaying  ’  file_name(i, : )  ’  ’  prep  ]) 
pause(l) 

if  size(file_name,l)  >  3 

rect  =  [(i-l)*105+5,5,200,200] ; 
else 

rect  =  C(i-l)*2i0+5,5,200,200] ; 

end 

hand  =  figureC ’Position’ .rect,  ... 

'lane' , [f ile_name(i, :)  ’  ’  prep  ], 

’ lumberTitle ’ , ’ of f  ’ ) ; 
if  i  ==  1 

display(aspectl .hand, auto) 
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elseif  i  =*  2 

display (aap«ct2 .hand, auto) 
elseif  i  ==  3 

display (aspact3, hand, auto) 
elsaif  i  ==  4 

display (aspsct4, hand, auto) 
elseif  i  ==  6 

display (aspsct 6 , hand, auto) 
elseif  i  ==  6 

display (aspsct 6, hand, auto) 
elseif  i  ==  7 

display (aspect7, hand, auto) 
and 

if  i  ==  1  I  datact  >=  1 
X  Indicata  targets 

for  j=l: langth(targats) 
if  targats(j.l)  "=  0 

text (targets (j , 1) , targets (j ,2) , 'X' ) 
and 
end 

and 

and 

and 

sat (ness , ' String ' , ' Output  Complete  * ) 
pause(B) 
figura(l) 
da 

hold  off 

X***************************************************************************** 
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%***************************************************************************** 
%  ROUT  HE  display  2.  b 

%*****************************************************♦*********************** 
Purpose:  Makes  encapsulated  postcript  file  of  all  displayed  figures  except 
for  figure  1,  the  nenu  figure. 

Called  by:  ear. a 

Quirks : 

1)  setup  for  aaxiaua  of  21  displayed  iaages  (22  figures) . 

%***************************************************************************** 

X  setup  aessage  window 
figure(l) 
cla 

hold  on 

axis( 'xy ’ , ’off’) 
mess  =  text(.l,.6,’  ’); 

X  check  to  see  that  figures  are  displayed 

h  =  figure( ’visible’ , 'off ’); 
if  h  <=  2 

set (aess, ’String’ , ’Figures  are  not  displayed  at  this  tine') 
pause (5) 
dell 

else 

delete (h) 
h  =  h  -  1; 
for  i=2:h 
figured) 

set  (aess, ’String’ , [’Printing  Figure’  i]) 
pause(l) 
if  i  ==  2 

print  fig2  -deps 
elseif  i  ==  3 

print  fig3  -deps 
elseif  i  ==  4 

print  fig4  -deps 
elseif  i  ==  5 

print  figB  -deps 
elseif  i  ==  6 

print  fig6  -deps 
elseif  i  ==  7 

print  fig7  -deps 
elseif  i  ==  8 

print  fig8  -deps 
elseif  i  ==  9 

print  fig9  -deps 
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alsail  i  ==  10 

print  liglO  -daps 
alsail  i  ==  11 

print  1  igll  -daps 
alsail  i  ==  12 

print  fig  12  -daps 
alsail  i  ==  13 

print  figl3  -daps 
alsail  i  ==  14 

print  ligl4  -daps 
alsail  i  ==  15 

print  liglS  -daps 
alsail  i  ==  16 

print  liglfi  -daps 
alsail  i  ==  17 

print  liglT  -daps 
alsail  i  ==  18 

print  ligl8  -daps 
alsail  i  ==  19 

print  ligl9  -daps 
alsail  i  ==  20 

print  lig20  -daps 
alsail  i  ==  21 

print  lig21  -daps 

end 

alsail  i  =*  22 

print  lig22  -daps 
and 
and 

and 

X  Finish  np 

sat (ness , ' String ' , * Output  Complete 1 ) 
pause (5) 
figured  ) 
cla 

hold  oil 
clear 


%***************************************************************************** 
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%************************************************************•***************•* 
X  ROUTIIE  dell.a 

X*****************»*************** ********************************************* 
X  Purpose :  Delates  the  current  set  of  figures  except  for  figure  1  which  is 
X  the  menu  figure  window. 

X 

X  Called  by:  sar.a,  displayl.a,  del2.a 
X 

y,***************************************************************************** 

h  =  figure( ’Visible’ , 'off '); 

for  i  =  h:-l:2 
delete(i) 
end 

clear  h  i 

figure(l) 

cla 

hold  on 

axis(’xy* , ’off’) 

X****************************************************************************** 
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X****************************************************************************** 
X  FUfCTIOI  display. a 

X****************************************************************************** 
Purpose:  Log  encodes  the  inage,  scales,  and  adjusts  the  intensity  of  the 
image  for  display  purposes. 

Called  by:  displayl.n 

Quirks : 

1)  Intensity  adjustments  are  hard-coded. 

■OTE: 

MatLab  displays  the  images  as  transposes,  z  said  y  coordinates  are 
reversed.  In  order  to  correct  this  so  that  target  coordinates  could 
be  plotted  over  the  image  and  natch,  the  image  is  transposed  just  before 
displaying . 

X***************************************************************************** 

function  display(X, fig, autoscale) 

X  Converting  to  <1B  and  scaling 
X  =  20*logl0(abs(X)) ; 

X  Selecting  maximum  and  minimum 

if  autoscale  <=  0  X  images  are  clipped  at  -SO  and  20  dB 
minimum  =  -50; 
maximum  =  20; 

elseif  autoscale  ==  1  X  images  are  ranged  between  mu  said  min  pixels 
minimum  =  min(min(X)); 
maximum  =  max (max (X) ) ; 
end 

X  scaling 

range  =  abs (maximum  -  minimum); 
if  range  *=  0 

X  =  round(((X  -  minimum) /rsuige) *255) ; 

else 

X  =  0  *  ones(size(X, 1) ,size(X,2)) ; 
end 

X  clipping  out  of  range  values 
if  autoscale  ==  0 
[i,j]  =  f ind(X<0) ; 
for  k  =  1  :  length(i) 

X(i(k) , j(k))  =  0; 
end 

[i.j]  =  f ind(X>255) ; 
for  k  =  1  :  length(i) 

X(i(k) , j(k))  =  255; 

end 
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•ad 


X  correcting  lor  Mat 1 aba  inaging  axis  transposition 

X  -  X'; 

X  Display  inage 
figure (lig) 

set dig, ‘PaperPosition’ , [1,3,6,5.75] , 'lextPlot ’ , 'new') 
hold  on 

axisC'xy', *ott*) 
colornap(gray) 
if  autoscale  ==  0 
inaged-75) 

elseif  autoscale  «»  -1 
inage(X-150) 
elseif  autoscale  =*  1 
iaage(X-llS) 
end 

hold  off 

%************************************eeeeee*e********************************** 
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%*******«********************************************************************** 
X  ROUTIIE  del2.a 

X****************************************************************e*******ee**** 
X  Purpose:  Oelstes  all  existing  .sat  and  .tax  files  in  current  directory 
X  that  would  be  created  by  a  session  of  ear. a 

X 

X  Called  by:  sar.u 
X  Calls:  dell .a 
X 

X***************************************************************************** 

if  exist (’unrotatedl. rat')  ==  2 
delete( 'unrotatedl .rat ’ ) 
end 

if  exist ('unrotated2. rat*)  ==  2 
delete ( 'unrotated2.rat ' ) 
end 

if  exist ('unrotated3. rat')  ==  2 
delete ( ’unrotated3 .rat ' ) 
end 

if  exi8t( ’unrotated4.rat')  ==  2 
delete( ’unrotated4 .rat ' ) 
end 

if  exist ('unrotatedS. rat')  ==  2 
delete ( 'unrotatedS .rat ’ ) 
end 

if  exist ( ' unrotat ed6 . rat ' )  ==  2 
delete( ’ unrot ated6. rat ' ) 
end 

if  exist('unrotated7.ut’)  ==  2 
delete ( ’unrotated7 .rat ' ) 
end 

if  exi8t( 'rotatedl .rat*)  ==  2 
delete ( 'rotatedl .rat ' ) 
end 

if  exist ('rotated2. rat')  ==  2 
delete ( *rotated2. rat') 
end 

if  exist ('rotated3. rat')  ==  2 
delete( 'rotated3.rat' ) 
end 

if  exist ( 'rotated4.rat')  ==  2 
delete ( 'rotated4 .rat ' ) 
end 

if  exist(’rotated5.rat')  ==  2 
delete ( 'rotatedS .rat  * ) 
end 

if  exist ('rotatedS. rat’)  ==  2 
delete ( 'rotatedS .rat ’ ) 
end 
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if  exist ('rotated7. nt')  ==  2 
delete (’rotated7. mat’) 

•ad 

if  exist ('rotatedd.aat')  ==  2 
delete ( 'rotatedd.aat') 
end 

if  exist ('detl.aat’)  ==  2 
delete( 'detl.aat  * ) 
end 

if  exist  ( '  det2 .  aat  ’ )  ==  2 
delete ( ' det2 . mat ' ) 
end 

if  exist (*det3. aat')  ==  2 
delete ( 'det3.aat ’ ) 
end 

if  exist ('averaged. aat')  ==  2 
delet e ( ' averaged . aat ' ) 
end 

if  exist ( 'spaaed. aat')  ■■  2 
delete ( ’ spaned .aat ’ ) 

•ad 

if  exist ('optiaal. aat')  ==  2 
delet • ( ' optiaal . aat ' ) 
end 

if  exist ( ' vhitened.aat' )  ==  2 
delete (' whitened . aat ’ ) 
end 

if  exist( 'statistics. aat')  ==  2 
delete ( ’ statistics .aat ' ) 
end 

if  exist ('statl.tex')  ==  2 
delete ( ’ statl . tex' ) 
end 

if  exist ('stat2.tex')  ==  2 
delete( ’ stat2 . tex ' ) 
end 

if  exist ('detl.eps')  ==  2 
delete ( *detl . eps ' ) 
end 

if  exist( 'det2.eps')  ==  2 
delete( 'det2.eps ’ ) 
end 

if  exist ('det3.eps')  ==  2 
delete( 'det3 . eps ' ) 
end 

if  exist (’files. aat')  ==  2 
delete ( 'files .aat') 
end 
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dell 


clear 

%•**•************************************************************************** 
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X****************************************************************************** 
X  &OUTIIE  show.stat.m 

X***************************************************************«************** 
X  Purpose:  Displays  the  statistical  characteristics  of  the  current  data  set 
X  If  statistics. mat  does  not  exist  sill  create  it.  Creates  two  .tex 

X  files  for  input  into  a  LaTeX  document ,  statl.tex  contains  the 

X  standard  deviation  to  aean  ratios  and  stat2.tex  contains  the 

X  target  to  clutter  ratios. 

X 

X  Called  by:  sar.a 
X  Calls :  stat . a 
X 

X***************************************************************************** 

if  exist (’ stat istics.aat’)  ~-2 
stat 

else 

load  statistics.aat 
load  files 
end 

X  setup  figure  window 
figure(l) 
cla 

hold  on 

axisCxy' ,  ’off) 

X  stain  title 

text(.2,l- 05, ‘Standard  Deviation  to  Mean  Ratio  (dB)’) 

subtitles  =  ['unrot  HH'  ;  ’unrot  HV*  ;  ’unrot  W’  ;  ... 

’rot  HH  ’  ;  ’rot  HV  ’  ;  ’rot  W  '  ;  . . . 

’averaged'  ;  'spaned  ’  ;  'optimal  ’  ;  'whitened'); 

X  sub  titles  and  output  of  standard  deviation  to  mean 
for  i  =  1  :  size(file_naae,l) 

text(  (i*0.15)-.l,l.Iil«_naae(i,:)  ) 
end 

for  i  =  1:10 

text(~. 15, l-( .05*i) ,subtitles(i, : )) 
end 

for  i  =  l:size(file_name,l) 
for  j  =  1:10 

text(  d*0.15)-.l,l-(-05*j),  nua2str(sdm(j ,i))  ) 
end 
end 

X  title 

text( .2, .45, 'Target  to  Clutter  Ratio  (dB)’) 

X  sub  titles  and  output  of  target  to  clutter 
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for  i  =  1  :  sized ile.nane,  1) 

text(  (i*0.15)-.l, .4,file_naae(i, :)  ) 
end 

for  i  =  1:10 

toxt(-.15, .4-(.05*i) ,subtitles(i, :)) 
end 

for  i  =  l:size(file_nane,l) 
for  j  =  1:10 

toxt(  (i*0. 15)-. 1, .4-(.06*j) ,  nu»2str(tc(j ,i))  ) 
end 
end 

X  print  to  ascii  filo  for  reading  by  LaTex 
fid  =  fopenCstatl.tex’  ,'n’); 
for  i  =  1:10 

fprintf (fid, ’Xs’ , subtitles(i, : )) ; 
for  j  =  l:size(file_naae, 1) 

fprintf (fid, '\t  ft  %4.2f 1 ,sda(i, j)) ; 
end 

fprintf (fid, ’\t  \\\\  Whline  \n’); 
end 

f  closed  id) ; 

fid  =  fopen(’ stat2.tex’ , ’»') ; 
for  i  =  1:10 

fprintf (fid. 'Xs' ,subtitles(i, :)); 
for  j  =  l:size(file_nane,l) 

fprintf (fid, *\t  *  X4.2f ' ,tc(i, j)) ; 
end 

fprintf  (fid,  *\t  \\\\  Whline  \n’ ) ; 
end 

fdose(f  id) ; 
clear 

%****************************************************************************** 
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%****************************************************************************** 
X  ROUTIIE  stat.B 

%***********#*****************************************************»*****«****** 
X  Purpose:  Creates  statistics. mat  which  contains  the  statistical  data  on  the 
X  current  image  set. 

X 

X  Called  by:  show.stat.a 
X  Calls:  targclut.a 

X 

X  Quirks : 

X  1)  Bill  not  give  information  on  files  that  have  not  been  created. 

X 

%***************************************************************************** 

time  =  fix (clock); 

X  setup  aessages 
figure(l) 
cla 

hold  on 

axis(’xy’.’off') 
aessl  =  text( . 1, .6, ’  ’); 
aess2  =  text (.1,. 5,'  ’); 
taess  =  text(.l,l, '  ’); 
aaess  =  text( .1, .95, '  ’); 

set (taess, 'String* , [’Starting  Time:  ’  nua2str(tiae(l,4))  ... 

nua28tr(tiae(l,5))  nua2str(tiae(l,6))  ]); 
set  (aaess, 'String' , 'Max  Run  Time:  17  ainutes');  pause(l) 


X  sda  =  standard  deviation  to  mean  ratio  (in  db) 

X  tc  =  target  to  clutter  ratio  (in  db) 

X 

X  matrices  have  the  following  fora 

X 

unrotated  HH  minus 
unrotated  HV  minus 
unrotated  VV  minus 
rotated  HB  minus 
rotated  HV  minus 
rotated  W  minus 
averaged  minus 
spaned  minus 
optimal  minus 
pwf  minus 

load  files 
X  unrotated 

if  exist ('unrotatedl. mat’)  ==  2 


X  unrotated  HH  plus 

X  unrotated  HV  plus 

X  unrotated  W  plus 

X  rotated  HH  plus 

X  rotated  HV  plus 

X  rotated  W  plus 

X  averaged  plus 

X  spaned  plus 

X  optimal  plus 

X  pwf  plus 


unrotated  HH  zero 
unrotated  HV  zero 
unrotated  W  zero 
rotated  HH  zero 
rotated  HV  zero 
rotated  W  zero 
averaged  zero 
spaned  zero 
optimal  zero 
pwf  zero 
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for  i  =  l:size(file_na*e,l) 

, 'String' , [’Loading  unrotated  ’  f ile_na»e(i, : )] ) 
paused ) 
if  i  ==  1 
load  unrotatedl 
elseif  i  ==  2 
load  unrotated2 
•laeif  i  ==  3 
load  unrotated3 
alaaif  i  ==  4 
load  unrotat «d4 
alaaif  i  ==  6 
load  unrotat edS 
alaaif  i  ==  6 
load  unrotatedd 
alaaif  i  7 
load  unrotatad7 
and 

aat(naaa2,  ’String’ Standard  deviation:  naan  ’  fila_nana(i, : )]  ,  ... 

’Visible’ , ’on’) 

panaa(l) 

adn(l,i)  =  20*logl0(abs(  std(std(hh))/nean(nean(hh))  )); 
ada(2,i)  =  20*logl0(abs(  std(std(hv))/nean(nean(hv))  )); 
sdn(3,i)  =  20*logl0(abs(  std(std(vv))/nean(nean(vv))  )); 

set(aasa2, ’String’ , [’ Target : clutter  ’  file_naae(i, :)] , 'Visible’ , ’on’ 
pause(l) 

tc(l,i)  =  targdut  (hh,  targets) ; 
tc(2,i)  =  targdut (hv, targets)  ; 
tc(3,i)  =  targdut  (vr,  targets  ) ; 
set(ness2, ’Visible’ , 'off’) 
pause(l) 


end 

end 

X  rotated 

if  exist (’rotatedl. nat’)  ==  2 
for  i  =  l:size(file_najne,l) 

set(nessl,  'String' ,  [’Loading  rotated  ’  file_name(i, : )] ) 
pause(l) 
if  i  ==  1 

load  rotatedl 
elseif  i  ==  2 
load  rotated2 
elseif  i  ==  3 
load  rotated3 
elseif  i  ==  4 
load  rotated4 
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elseif  i  ==  5 
load  rotatedS 
alsaif  i  ==  8 
load  rotatedS 
alsaif  i  ==  7 
load  rotatad7 
and 

sat (mess2, ’String* , ['Standard  deviation: naan  ’  f ile_name(i, : )] , _ 

’Visible’ , 'on') 

paused  ) 

sdm(4,i)  =  20*logl0(abs(  std(std(hh))/mean(mean(hh))  )); 
sdm(5,i)  =  20*logl0(abs(  std(std(hv))/mean(mean(hv))  )); 
sdm(6,i)  =  20*logl0(abs(  std(std(w))/mean(aeaa(vv))  )); 

set(mess2, ’String’ , [’Target: clutter  ’  file_nana(i, :)] , ’Visible’ , ’on’) 
pause(l) 

tc(4,i)  =  targclut (hh, targets ) ; 
tc(5,i)  =  targclut (hr, targets ) ; 
tc(6,i)  ■  targclut (vv, targets ) ; 
set(mess2, ’Visible’ , ’off’) 
pause(l) 


and 

and 

dear  hh  hv  vr 
X  averaged 

if  exist (’averaged. mat’)  ==  2 

set (mess 1, ’String’ , 'Loading  averaged  ’) 
pausa(l) 
load  averaged 

for  i  =  1  :  size(file_name,l) 
if  i  ==  1 
temp  =  aspectl; 
alsaif  i  ==  2 
tamp  =  aspect 2; 
alsaif  i  ==  3 
temp  =  aspect 3; 
elseif  i  ==  4 
temp  =  aspect4; 
alsaif  i  ==  5 
tamp  =  aspects; 
alsaif  i  ==  8 
temp  =  aspect6; 
alsaif  i  ==  7 
temp  =  aspect7; 
and 

sat  (mass2l  'String’ ,  ['Standard  deviation  mean  '  file_name(i, :)]  ,  ... 
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'Visible' , ’on’) 


paused) 

sda(7,i)  »  20*logl0(abs(  std(std(teap))/aean (Bean (teap))  )); 
set(aess2, 'String' , ['Target: clutter  *  lile_naae(i, : )] , 'Visible' , 'on') 
panse(l) 

tc(7,i)  =  targclut(teap, targets); 
end 

set(mess2, 'Visible' .'oil') 
panse(l) 

end 

X  span 

if  exist('spaned.aat')  ==  2 
set(aesal,  'String' ,  ’Loading  spaned  ') 
paused  ) 
load  spaned 

lor  i  =  1  :  size(lile_naae,l) 
if  i  ==  1 

teap  =  aspect 1 ; 
elseif  i  ==  2 
teap  =  aspect2; 
elseil  i  ==  3 
teap  =  aspect 3; 
elseil  i  ==  4 
teap  =  aspect4; 
elseil  i  ==  5 
teap  =  aspects; 
elseil  i  ==  6 
teap  =  aspects; 
elseil  i  ==  7 
teap  =  aspect7; 
end 

set(aess2, 'String' , ['Standard  deviation: mean  '  lile_name(i, : )] ,  ... 

'Visible' , 'on') 

pause(l) 

sdn(8,i)  =  20*logl0(abs(  std(std(teap))/aean(aean(t«Bp))  )); 
set(aess2, 'String' , ['Target: clutter  ’  file_name(i, : )] , 'Visible' , 'on') 
paused) 

tc(8,i)  =  targclut(teap, targets); 
end 

8et(aes82, 'Visible' , 'oil') 
pause(l) 

end 

%  optiaal 

il  exist (' opt iaal.aat')  ==  2 

set (aessl , ’ String ' , 'Loading  optiaal  ' ) 
pause(l) 
load  optiaal 
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for  i  =  1  :  a ize (f ile.naae , 1 ) 

if  i  ss  i 

tenp  =  aapactl; 
alsaif  i  ==  2 
teap  s  &apect2; 
alsaif  i«3 
taap  =  aspact3; 
alsaif  i  ==  4 
taap  =  aspact4; 
alsaif  i  =®  5 
taap  ®  aspects; 
alsaif  i  ==  8 
taap  =  aapsctfi; 
alsaif  i  ==  7 
taap  ®  aspect 7; 
and 

set(aeas2, ’String’,  ['Standard  deviation :aaan  ’  f  ile_naae(i, : )]  ,  ... 

'Visible', 'on') 

pause(l) 

sda(9,i)  =  20*logl0(aba(  atd(std(teap))/aean(nean(taap))  )); 
8at(aaaa2, 'String' , ['Target : clutter  '  file_naae(i, : )] , 'Visible' ,  'on') 
pausa(l) 

tc(9,i)  ®  targclut (taap, targets); 
end 

aat(aa882, 'Visible' , 'off') 
pauseCl) 

end 
X  pwf 

if  eziat('sbitened.Bat')  ==  2 

setCaesai, 'String', 'Loading  whitened  ') 
paused) 
load  whitened 

for  i  =  1  :  size(file_naae,l) 
if  i  «  1 
teap  =  aapactl; 
elaeif  i  ==  2 
taap  =  aspect2; 
elaeif  i  ==  3 
teap  =  aspect3; 
elaeif  i  ==  4 
teap  =  aspect 4; 
elaeif  i  ®»  5 
teap  =  aspects; 
elaeif  i  ==  6 
teap  =  aapect6; 
elaeif  i  ==  7 
teap  =  aspect7; 
end 

set(aeas2,  'String',  [’Standard  deviation:aean  ’  file_naae(i, :)]  ,  ... 
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’Visible’ , 'on*) 

paused  ) 

sda(10,i)  =  20*logl0(abs(  std(std(teap))/*ean(uean(te*p))  )); 
set (aess2 , 1 String Target : clutter  ’  f ile.naae ( i , : )] , * Visible ’ , ’ on ’ ) 
pause(i) 

tc(10,i)  =  targdut(teap,  targets); 
end 

set(aess2,  'Visible* ,  ’off) 
pause(l) 

end 

clear  aspectl  aspect2  &spect3  aspect 4  aspects  aspects  aspect? 

%  wrap  up 

set (aessl, ’String* , 'Saving  as  statistics .aat * ) 
pause(l) 

save  statistics  sda  tc 

set  (aessl,  ’String’ ,  ’Statistical  Routine  Coaplete’) 
pause(S) 
figured) 

da 

hold  off 

X*********************************************************************«******** 
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x****************************************************************************** 
X  FOICTIOI  targclut.a 

X****************************************************************************** 
X  Purpose:  calculates  the  target  to  clutter  ratio  of  an  image,  given  the 
X  known  coordinates  of  targets 

X 

X  Called  by:  stat.m 
X 

X  Quirks : 

X  1)  assumes  a  "large"  image  with  small  number  of  targets. 

X 

X***************************************************************************** 

function  Cy]  =  targdut(x, targe) 

X  z  -  input  image 
X  targe  -  target  coordinates 

r  =  5;  X  5  chosen  to  get  good  average  over  target,  but  not  include  clutter. 

X  lots  that  in  the  following  algorithm,  if  target  is  in  corner  it  sill  not 
X  be  used 
targr  =  □; 

for  i  =  1 : size (targe, 1) 
temp  =  □ ; 

if  targc(i,l)  '=  0  ft  targc(i,l)>=r+l  ft  targc(i,2)>=r+l  ft  ... 
targc(i,l)<=8ize(x,l)-r-l  ft  targc(i,2)<=size(x,2)-r-l 
temp  =  x(  targc(i,l)-r:targc(i,l)+r  ,  targc(i,2)-r:targc(i,2)+r  ); 
elseif  targc(i.l)  '=  0  ft  targc(i,2)>=r+l  ft  ... 

targc(i,l)<=size(x,l)-r-l  ft  targc(i,2)<=size(x,2)-r-l 
temp  =  x(  l:targc(i,l)+r  ,  targc(i,2)-r:targc(i,2)+r  ); 
elseif  targc(i.l)  '=  0  ft  targc(i,l)>=r+l  ft  ... 

targc(i,l)<=size(x,l)-r-l  ft  targc(i,2)<=size(x,2)-r-l 
temp  *  x(  targc(i,i)-r:targc(i,l)+r  ,  l:targc(i,2)+r  ); 
elseif  targc(i,l)  '=  0  ft  targc(i,l)>=r+l  ft  targc(i,2)>=r+i  ft  ... 
targc(i,2)<=size(x,2)-r-l 

temp  =  x(  targc(i,l)-r:size(x,l)  ,  targc(i,2)-r:targc(i,2)+r  ); 
elseif  targc(i.l)  ”*  0  ft  targc(i, l)>=r+l  ft  targc(i,2)>=r+l  ft  ... 
targc(i, l)<=size(x, l)-r-l 

trap  =  x(  targc(i, l)-r:targc(i, l)+r  ,  targe (i, 2) -r: size (x, 2)  ); 
end 

targr  =  [targr  ;  reshapeCtemp,  size(temp,l)*size(temp,2)  ,1)  ]; 
end 

y  =  20*logl0(abs(  mean(mean(targr))/mean(mean(x))  )); 

X  mean  for  clutter  was  taken  over  entire  image  to  simplify  routine 
X  contribution  of  target  to  entire  image  should  be  small  if 
X  number  of  targets  is  small  and  image  is  large 

X****************************************************************************** 
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%****************************************************************************** 

X  aotrriiB  deii.a 

X****************************************************************************** 
X  Purpose:  Oalatea  the  currant  aat  of  figures  axcapt  lor  figure  1  which  is 
X  tha  menu  figure  window. 

X 

X  Called  by:  aar.n,  display 1.*,  del2.a 
X 

X***************************************************************************** 

h  *  fignre('Visible’, ’oil’); 

lor  i  =  h:-l:2 
delete (i) 
and 

dear  h  i 

ligure(l) 

da 

hold  on 

axis('xy' , 'oil') 

X****************************************************************************** 
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X****************************************************************************** 
X  FUICTIOI  detection. a 

X*****o**********^******e*****************e************************************ 
X  Purpose:  Selects  the  fora  of  the  detection  and  displays  coordinates  of 
X  known  and  detected  targets. 

X 

X  Called  by:  sar.a 

X  Calls:  detectl.a,  detect2.n 

X 

X  Quirks : 

X  1)  lot  set  up  to  handle  unrotated  iaages. 

X  2)  Cut-off  value  for  determining  threshold  value  is  hard-coded. 

X 

X***************************************************************************** 

function  □  =  detect ion (type) 

tiae  =  fix(dock); 

X  setup  aessage 
figure(l) 
hold  on 

&xis(’xy* , 'off') 
cla 

ness  «  text(.l, .6, ’  '); 
taess  *  text(.l,l,'  '); 
aaess  a  text( . i, .95. '  »); 

set  (taess,  'String*,  [’Starting  Tiae:  *  nua2str(tiae(l,4))  ... 

nua2str( tiae (1,5))  »:»  nua2str(tiae(l,6))  ]); 
set (aaess, ’String’ , 'Max  Run  Tiae:  indeterainant ’ ) ;  pause(l) 

load  files 

if  typo  ==  0 
detect  =  0; 

set (aess, ’String’ , 'lo  Detection  Done’) 
pause(5) 
figured) 
cla 

hold  off 
elseif  type  ==  i 
detect  a  i; 
type  =  -1; 
elseif  type  as  2 
detect  =  2; 
type  =  -1; 
elseif  type  ==  3 
detect  =  3; 
type  =  -1; 
end 
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ufi  files  director;  file.nane  prop  detect 

if  type  <  0 

if  prep  *=  'rot* 

set  (seas , 'String' , 'Loading  rotated' ) 
paase(l) 

if  exist ('rotatedd.aat')  ==  2 
load  rotatedd.aat 

else 

if  exist ('rotatedl.aat')  ==  2 
load  rotatedl 
aspect  1  =  hit; 

end 

if  exist ('rotated2. aat’)  ==  2 
load  rotated2 
aspect 2  =  hh; 
end 

if  exist('rotated3.aat*)  ==  2 
load  rotated3 
aspect3  a  hh; 
end 

if  exist ('rotatedd.aat')  --  2 
load  rotatedd 
aspeetd  a  hh; 
end 

if  exist (’rotatedS. aat’)  as  2 
load  rotatedS 
aspects  a  hh; 
end 

if  exist (’rotatedS. aat’)  ==  2 
load  rotatedS 
aspects  =  hh; 
end 

if  exist('rotated7.aat')  ==  2 
load  rotated7 
aspect 7  a  hh; 

end 

save  rotatedd  aspect 1  aspect2  aspect3  . . . 

aspeetd  aspects  aspects  aspect7  targets 

end 

elseif  prep  ==  'are' 

set  (boss,  'String' ,  'Loading  averaged') 
panse(l) 
load  averaged 
elseif  prep  ==  'spn' 

set  (boss  ,  'String' ,  'Loading  spaned') 
panse(l) 
load  spaned 
elseif  prep  ==  *o«s’ 

set  (boss, 'String' , 'Loading  optiaal’) 
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paused) 
load  optimal 
elseif  prep  ==  ‘pvt* 
set  (mess,  ’String’ ,’  Loading  whitened') 
panse(l) 
load  whitened 
end 

if  detect  ==  1 

set  (mess,  'String' ,  [’Multiple  Aspect  Detection  on  ’  prep  '  (scaled  threshold)']) 
panse(l) 
detects 

elseif  detect  ==  2 

set (aess, 'String' , [’Multiple  Aspect  Detection  on  '  prep  '  (preset  threshold)']) 
pause(l) 
detect2 

elseif  detect  ■■  3 

set (mess, 'String* , [’Subtraction  on  ’  prep]) 
pause(l) 
detect2 
end 

X  Output  information 
fignre(l) 
cla 

if  prep  ==  'non' 

text (0,1. 05, 'Must  select  a  pre-processed  form!’) 

pause (5) 

cla 

end 

end 

%****************************************************************************** 


B-45 


%****************************************************************************** 
X  FU1CTI0I  detect 1 .■ 

X****************e************************************************************* 
X  Purpose:  Subtracts  two  aspects  froa  each  other. 

X 

X  Called  bj:  detect ion. a 
X  Calls:  detect2.a 
X 

X***************************************************************************** 

function  [y]  =  detectl (aspect 1 , aspect2) 

X  difference  between  aspect  1  and  aspect2 
y  =  (  abs(aspectl)  -  abs(aspect2)  ); 

X  adding  minianm  value  back  in 
ainiaua  =  0.001; 
y  *  y  +  min(ain(y))  +  ainiaua; 

X****************************************************************************** 
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%******************************************************************«*********** 
X  ROUT I IE  detect2.m 

X****************************************************************************** 
X  Purpose:  Uses  threshold  detection  on  multiple  images. 

X 

X  Called  by:  detection. m 
X 

X  Calls:  thresh. ■ 

X 

X  Quirks : 

X  1)  Radius  for  grouping  pixels  for  a  single  target  is  hard-coded. 

X  2)  Cutoff  values  are  hard-coded. 

X 

X***************************************************************************** 


upper  *  1.0 
inc  =  .05 
lower  =  .05 


X  upper  limit  of  cutoff 
X  incrementing  cutoff 
X  lower  limit  of  cutoff 


X  subtraction  detection 
if  detect  ==  3 

set (mess , ’ String ’ , Subtraction  on  '  prep] ) 
pause(l) 

if  size(file_name,l)  >=  1 
aspect 1  =  aspect 1; 
end 

if  size(file_name,l)  >=  2 

aspect2  =  det act 1 (aspect 1, aspect 2) ; 
end 

if  size(file_name,l)  >=  3 

aspect3  =  detect 1 (aspectl, aspect3) ; 
end 

if  size(file_name,l)  >=  4 

aspect4  =  det act 1 (aspect 1, aspect4) ; 
end 

if  size(file_name,l)  >=  5 

aspects  =  detectl(aspectl, aspects); 
end 

if  size(file_name,l)  >=  6 
aspects  =  det ectl( aspectl .aspects) ; 
end 

if  size(file_name.l)  >=  7 

aspect 7  =  detect 1 (aspectl .aspect 7) ; 
end 

set  (mess,  ’String’ ,  ['Subtraction  Detection  on  '  prep]) 
pause (1) 

end 


detected  =  □ ; 

for  cutoff  =  lower: inc: upper 

X  set  threshold  using  most  positive  aspect 
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if  detect  ==  2 

threshold  =  (  max (max (abs( aspect 1)))  +  min(min(abs (aspect 1)))  )*cutoff; 
end 

X  Ron  detection  on  all  aspects 
if  size(f ile_naae, 1)  >=  1 
if  detect  ==  1  I  detect  ==  3 

threshold  =  (  irax(max(abs(aspectl)))  +  Bin(min(abs(aspectl)))  )*cutoff; 
end 

targsl  =  thresh(aspectl .threshold) ; 
if  detect  '=  3 
targsc  =  targsl; 
end 
end 

if  size(f ilejaane, 1)  >=  2 

if  detect  ==  1  |  detect  ==  3 

threshold  *  (  aax(Bax(abs(aspect2)))  +  min(ain(abs(aspect2)))  )*cutoff; 
end 

if  detect  "=  3 

targsc  =  [targsc  ;  thresh (aspect 2, threshold)] ; 

else 

targsc  =  thre8h(aspect2, threshold); 
end 
end 

if  size(file_name,l)  >=  3 

if  detect  ==  1  |  detect  ==  3 

threshold  =  (  max (max (aba (aspect 3)))  +  min(min(abB (aspects)))  )*cutoff ; 
end 

targsc  =  [targsc  ;  thresh(aspect3 , threshold)] ; 
end 

if  size(file_name,l)  >=  4 

if  detect  ==  1  I  detect  ==  3 

threshold  =  (  max(max(abs(aspect4)))  +  min(min(abs(aspect4)))  )*cutoff; 
end 

targsc  =  [targsc  ;  thresh(aspect4, threshold)] ; 
end 

if  sizetfile.name.l)  >=  S 

if  detect  ==  1  |  detect  ==  3 

threshold  =  (  max(max(abs(aspect5)))  +  min(min(abs (aspects)))  )*cutoff; 

end 

targsc  =  [targsc  ;  thresh(aspect5, threshold)] ; 
end 

if  size(file_naae,l)  >=  6 

if  detect  ==  1  I  detect  ==  3 

threshold  =  (  max(aax(abs(aspect6)))  +  min(min(abs(aspect6)))  )*cutoff; 
end 

targsc  =  [targsc  ;  thresh(aspect6, threshold)] ; 
end 

if  size(file_naffle,l)  >=  7 

if  detect  ==  1  I  detect  ==  3 

threshold  =  (  maz(max(abs(aspect7)))  +  min(min(abs(aspect7)))  )* cut off ; 
end 
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targsc  =  [targsc  ;  thresh(aspect7 , threshold)] ; 
and 

X  For*  composite  image  with  targe  and  then  detect 

temp  *  0  *  ones(  size(aspectl,l)  ,  size(aspectl ,2)  ); 
lor  j  =  1: size (targsc, 1) 

temp(  targsc(j ,1)  ,  targsc(j,2)  )  =  255; 
end 

targsc  =  thresh(temp.lOO); 
clear  j  temp 

X  Determine  the  number  of  false  alarms  and  targets  detected  store  in  matrix 
lal  =  0;  lac  =0;  dl  =  0;  dc  =  0; 
r  =  16;  X  radius  around  target  coordinates 
lor  i  =  l:size(targsl,l) 
lor  j  =  l : size (targets, 1) 

il  targsl(i, : )<=targets(j , :)+r  *  targsl(i, : )>=targets(j , :)-r 
dl  =  dl  +  1; 
end 
end 
end 

lor  i  =  1 : size (targsc , 1) 
lor  j  =  l:size(targets,l) 

il  targsc(i, :)<=targets(j , :)+r  k  targsc(i, :)>=targets(j , :)-r 
dc  =  dc  +  1; 
end 
end 
end 

lal  a  8ize(targsl,l)  -  dl; 
lac  =  size (targsc, 1)  -  dc; 

detected  a  [  detected  ;  cutoll  size (targets, 1)  lal  dl  lac  dc] ; 

end 

X  print  to  ascii  lile  lor  possible  future  use 
lid  =  fopen(’det.tex’ , ’»’); 
for  i  =  1 : size (detected, 1) 

1 or  j  =  l:size(detected,2) 

IprintKfid,  ’X4.21  \t’  ,detected(i,  j))  ; 
end 

IprintKfid, ’\t  \n’); 

end 

lclose(lid) ; 
il  detect  ==  3 

save  det3  detected  targets  . . . 

aspect 1  aspect2  aspect3  aspect4  aspects  aspects  aspect7 

end 

X  plot 
dell 

il  detect  ==  1 

figured Position’ ,[5,5,515,400] .... 

’lame’ , [’Scaled  threshold  detection  on  ’  prep] , ’ lumberTitle ’ , ’oil’) ; 
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elseif  detect  ==  2 
figureC 'Position' ,[5,5,515,400] .... 

’■ame* , [ 'Preset  threshold  detection  on  ’  prep] , 'lumberTitle' , 'off ') 
elseif  detect  ==  3 

figure ('Poe it ion' , [5,5,515,400] .... 

' lane ’ , [ ’ Subtract ion  detect ion  on  ’  prep] , ’ lumberTitle ' , ’ of f ' ) ; 

end 

plot(detected(: ,l),detected(: ,3) , 'w-. ' .... 
detectedC : , 1) , detect ed( : ,4) , 'w — ' ,  ... 
detectedC: ,1) , detectedC : ,5) , 'u: ’ ,  ... 
detectedC: ,1) .detectedC: ,6), 'w-' ,  ... 
min(min(detected( : ,1))) , (max(max(detected))+l) , 'k. ' ,  ... 
detectedC: ,1) .detectedC: ,2)+0.1, ’r-’ ,  . . . 
detectedC: ,1) .detectedC: ,2)-0. 1, ’r-’) 

X  use  next  six  lines  to  set  axis  for  comparison  of  plots 
xmin  =  0; 
xnax  =  1.0; 
ymin  =  0; 
ymax  =  20; 

axisC [xnin , xnax , ymin , ymax] ) 
set .axis  =  1; 
xlabelC 'Threshold' ) 

ylabelC' lumber  of  Detections  or  False  Alarms') 

X  move  these  statements  to  just  before  the  clears 
X  in  order  to  print  the  legend  with  the  figure 
if  detect  ==  1 
print  detl  -depa 
elseif  detect  ==  2 
print  det2  -deps 
elseif  detect  ==  3 
print  det3  -deps 
end 


if  set.axis  ==  0 
xmin  =  minCdetectedC : ,1)); 
max  =  maxCdetectedC:  ,D); 

ymin  =  minCminCdetectedC : ,2:sizeCdetected,2)))); 
ymax  =  maxCmaxCdetectedC : ,2:size(detected,2)))); 
end 

x  =  Cxaax  +  xnin)  *  .5; 
dx  =  Cxaax  -  xmin)  *  .1; 
y  =  (ynax  +  ymin)  *  .75; 
dy  =  (ymax  -  ymin)  *  .05; 

legend  =  [  SAR  false  alarms  ’  ;  ... 

’ -  SAR  detected  ’  ;  ... 

' . VASAR  false  alarms  ’  ;  . . . 

’ _  VASAR  detected  ']; 

text(x,y, [num2str(size(targets, 1))  '  known  targets  ']) 
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lor  i  =  1: size (legend, 1) 
if  i  ==  1 

text(x,y+(dy*i)+(dy*.5) ,legend(i, 1:10)) 
elseif  i  ==  2 

text(x,y+(dy*i),l«gend(i,l:10)) 

else 

text (x , y+(dy*i)+(dy* . 25) , legend ( i , 1 : 10) ) 

•nd 

text (x+dx , y+ (dy*x) , legend ( i , 1 1 : size (legend , 2) ) ) 
end 

X  Move  print  stateaents  to  here  to  print  legend  with  figure 

clear  upper  loner  cutoff  threshold  targsl  targsc  fal  dl  fac  dc 
clear  xaax  xain  yn&x  yain  x  dx  y  dy  legend  ans  i 

%************♦***♦*****♦********♦********************************************* 
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X****************************************************************************** 
X  FUICTIOI  thresh. m 

Xee**************************************************************************** 
Purpose:  Threshold  detection. 

Called  by:  detect2.m 

Quirks: 

1)  Radius  for  grouping  pixels  for  a  single  target  is  hard-coded. 
X***************************************************************************** 
function  [targs]  =  thresh (aspect 1,  threshold) 

X  detection 

[x,y]  =  f ind(abs(aspectl)  >=  threshold); 
targs  =  [x  y]; 

X  group  returns  within  specified  radius  and  average 
r  =  15;  X  radius  of  pixels 
newtargs  =  □ ; 
for  i  =  1: size (targs ,1) 
temp  =  □ ; 

for  j  =  1 : size (targs , 1) 

if  targs(j,l)<=(targs(i,l)+r)  *  targs(j,l)>=(targs(i,l)-r)  ft  ... 
targs(j ,2)<=(targs(i,2)+r)  ft  targs(j,2)>=(targs(i,2)~r) 
temp  =  [  temp  ;  targs(j,:)  ]; 
end 
end 

temp  =  [round (mean (temp(: ,1)))  round(mean(temp( : ,2)))] ; 
newtargs  =  [  newtargs  ;  temp  /  length(j)]; 
end 

targs  =  newtargs; 

X  Get  rid  of  duplicate  coordinates 
for  i  =  1 : size (targs, 1) 
for  j  =  i : size (targs, 1) 
if  i  •=  j 

if  targs( j , l)<=(targs(i , l)+r)  ft  targs(j ,l)>=(targs(i,l)-r)  ft  ... 
targs(j ,2)<=(targs(i,2)+r)  ft  targs(j ,2)>=(targs(i,2)-r) 
targs (i, : )  =  [00]; 

end 

end 

end 

end 

for  i  =  size(targs,l) : —1 : 1 
if  targs (i, 1)==0 
if  i  *=  1 

targs  =  [  targs(l:i-l, :)  ;  targs (i+ 1 : size (targs, 1) ,: )  ]; 

else 

targs  =  targs (i+l:size(targs,l),:); 
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•ad 

•ad 

•ad 

%****************************************************************************** 
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